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Abstract 

This lecture addresses the concept of form factor in magnetic scattering of thermal 
neutrons, analyzing its meaning, discussing its measurement by polarized neutrons and 
detailing its computation for the ions by the spherical tensor operator formalism. 

1 INTRODUCTION 

It is an experimental fact that the magnetic scattering of thermal neutrons from a magnetic 
material can be detected only over a fraction of the accessible range of the scattering vector and, 
unlike the nuclear scattering, vanishes out as this gets large. Ascribed to the spatial distribution 
of the magnetic densities, which is modulated at the atomic scale and thus at the scale of the 
wavelength of the thermal neutrons, the underlying variation in the scattering amplitudes 
is extracted out through form factors associated to ions. A concrete advantage is that these 
depend on radial integrals that can be computed at once and tabulated for each ion and valence 
state. Within the so-called dipole approximation the only unknowns then are the proportions 
of the spin and orbital magnetic moments attached to the ions. A further interest is that the 
expression of the differential scattering cross section is considerably simplified and adapted for 
the quantitative analysis of the static and dynamic correlation functions of magnetic moments. 
A number of magnetic phenomena thus can be investigated with more ease, for instance the 
long range or short range magnetic orders with the associated magnetic transitions or crossovers 
and the magnon excitations emerging from these orders, the configurations and dynamics of 
magnetic moments in spin glasses or in spin fluids, the excitation spectrum of low dimensional 
magnets or of arrays of molecular nanomagnets, et cetera. 

With polarized neutrons, measurements can be performed at precisions which, interestingly, 
allow probing the actual spatial distribution of the magnetic densities. A wholly new set of 
features associated to these then become available to analysis, for instance the quantum states 
of the concerned ions, which requires the exact computation of the associated form factors, 
but also the degree of covalency with neighboring ions or of hybridization with conduction 
electrons or else the itinerant nature of the magnetic densities, which all might question about 
the use itself of form factors associated to single ions, the spin derealization and polarization 
effects, which might inform about the exchange paths and mechanisms in magnets, et cetera. 

We qualitatively examine in the following the concept of form factor, then describe the 
method used for its precise measurement by polarized neutron in collincar magnets and finally 
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provide with a survey of the algebra involved for their computation by the spherical tensor 
operator formalism. The aim of the lecture was to refresh our mind with concepts and method- 
ologies that are not new but unfamiliar because they are not of frequent use. Accordingly, there 
is nothing in this manuscript that cannot be found elsewhere. Nevertheless, no attempt will be 
made to exhaustively account for the literature, which instead will be limited to the minimum. 



2 SCATTERING AMPLITUDES 

The scattering amplitude A for a neutron-sample collision process, where the neutron spatial- 
spin incoming state |fcjfj) is transfered into the spatial-spin outgoing channel \kfVf) while the 
sample initial state |Aj) is transformed into the final state |Af), is computed equal to the matrix 
element 



A 



{(h^Xi) - {(fc>/)A/)}] - ((fc>/)A/) | V | (km)\i)) (2.1) 

in the Born approximation [1]. m n is the neutron massQ- V is the neutron-sample interaction 
potential. 

A first contribution to V is the neutron-nucleus nuclear interaction potential. Since its 
range (~ 1.5 x 10 -1 m) and the nucleus radius (~ 10 -14 m) are several orders of magnitude 
smaller than the wavelength of thermal neutrons, it essentially gives rise to isotropic s-wave 
neutron scattering and thus is well approximated by the Fermi-Dirac pseudo-potential 

V N = -Vb p i(rVR P ) (2.2) 



where S is the Dirac delta generalized function, r n the neutron position operator and R p the 
p-th nucleus position operator. b p = A p l t ®<7 + B p i p • er, for the nuclear interaction between 
two nucleons depends on their spin. 1 L ^ a is the unit operator in the p—th nucleous (g> neutron 
spin space, t p the p-th nucleus spin operator and i<x the neutron spin operator. The scattering 
lengths A p and B p are complex numbers that vary from nucleus to nucleus and as a function 
of the incident neutron energy, with the imaginary part representing radiative capture [2] . 

A second contribution to V, which naturally emerges from the non relativistic limit of the 
neutron Hamitonian in an electromagnetic field [1] , is the interaction potential 

V M = -fi n ■ B(r„) = -7 n /i N <x • B(r n ) (2.3) 

of the neutron magnetic moment p, n = 7„ /^n <t with the magnetic induction B(r n ) created by 
the sample at the neutron position r n . 7 n = —1.91348 is the neutron gyromagnetic ratio and 
/in the nuclear magneton q If the current density j(r) engendering the magnetic induction 
B(r n ) is stationary then 

B(?n) = ^ ! m A F "~! dr (2.4) 
4?r ./o r„-r P 



1 m n = 1.67492 10" 27 kg (939.57 MeV/c 2 ) = 1838.62 m c = 1.001375 m p , where m c = 0.910956 10" 30 kg 
(= 0.511 MeV/c 2 ) is the electron mass and m p = 1.67261 10~ 27 kg (= 938.28 MeV/c 2 ) the proton mass, c = 
2.997925 10 s ms -1 is the light speed and h = A = 1.05459 10~ 34 Js the Dirac constant (Planck constant/27r). 

2 m = ML = 5.05095 10" 27 JT" 1 . e (electron charge) = -1.60219 10 19 C. 
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where ^.q is the vacuum permeability d and Q any volume of the space with a boundary dtt 
strictly enclosing (cf. eq. (|A.7p and eq. (|A.8p ), An external magnetic induction also influ- 
ences the neutron spatial-spin state, which allows guiding or rotating the neutron polarization 
or else, in the case of strong gradients, channelling the neutron, but it affects the scattering 
by modifying the neutron statistical ensemble and not through the neutron-sample scattering 
amplitude. 

A neutron also reacts to the electrostatic field E(r n ) of the nuclear and electron charge 
densities through the first order relativistic corrections to the neutron Hamiltonian, namely 
the Spin - Orbit (SO) + Darwin Correction (DC) interaction potentials 

V so + V DC = a ■ (E(f n ) A p n ) - W ?n • E(r n ) (2.5) 

but these provide with negligibly small scattering amplitudes as compared to those nuclear 
(V N ) and magnetic (V M ) and can be safely ignored [2]. A little attention should be paid only 
to the N-SO elastic interference, bringing about, for initial neutron polarization Pj, a scattered 
intensity oc Pj • (kf A fcj) and, for null Pj, a scattered polarization ockj/\ fej, which might give 
rise to detectable effects when the samples shows small nuclear coherent scattering, contains 
strongly absorbing species or possesses a non centrosymmetric crystal structure [2]. 

The integration over the neutron space variable in the scattering amplitudes is straightfor- 
ward, since (u n \ki) — exp{i(/cj • u n )} and (fc/|s n ) = exp{— i(kf • s n )} are plane waves. Appropri- 
ately inserting the closure relation J \r n )(r n \ df n = J \s n )(s n \ ds n = J \ u n )(u n \ du n = • • • = 1 
and making use of the Kronecker symbol 5a,b to express that Sa,b = 1 A = B and 

8am = A 7^ _B, we get 

• for the nuclear scattering amplitude 

vn — < -> 

-^((kfVf)Xf)\VK\(k i u i )\ i )) = 

- ^2 {A p 6x f ,\i Sv f ,vi + Bp (A/| t R \Xi) ■ (uf\ <x |^)} exp{i(x-R p )} (2.6) 
p 

bearing in mind that (s n \ <5(r n — R p ) \u n ) — (s n \u n ) 5(u n — R p ) = S(s u — u n ) S(u n — R p ), 

• for the magnetic scattering amplitude 

7TL I — * — * \ —* 

- \( k f"f) x f) I v m I {km)Xi)J = {vf I it I Vi) ■ Sx f ,xt(ii) 

£\ f ,\i(£) = 7n r 7: 1 ? iZA(\f [ j(r) exp{i(x- f^} d? A;\ (2.7) 

2/l B ^ \ JQ I 

by using the identity Zg = -V?(~) = -V?(^s J exp{i(x • r)} dx) (cf. eq. (|A.2j) ). 

k = ki — kf defines the scattering vector, is the Bohr magneton and ro the classical radius 
of the electron 0. x and r are ordinary vectors so that 

ixA(A/| / J(f)exp{i(x • r)} dr |A<) = ix A / (A/| j(f) |A<) exp{i(x • r)} df (2.8) 



3 



fx Q = 4tt 10 



A m-1 

tt 5 - = -9.27410 10" 24 JT" 1 . r = = 2.81794 10" 
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The vector field (A/| j(r) |Aj) can always be split into the sum of an irrotational component 
ji \.(r) = — Vf (Xf \ S(f) \Xi) and a solenoidal component A .(r) = Vp A (A/| M(f) |Aj) : 



(A/| j(f) |Ai) = -V f (A/| H(r) |A*> + V ? A (A/| M(r) |A*> 



(2.9) 



where 



(A/| S(f) |A,> = - 



(A/ 1 M(f) | A,) = — 



V s ~ - (A/ 1 j(s) \Xi) 



ds 



V s ~ A(A/| |Ai) 



(2.10) 



(cf. eq. (|A.9|1 ). The irrotational component x^f) does not at all contribute to the scattering 
amplitude, since A j{ f x . (f) = so that —ikI\ J n ji ^ . (r) exp{i(i?-r)} dr = 0. On the con- 
trary, the solenoidal component A (r) is /uHy probed by the neutron, since Vf? • j£ A . (r) = 
so that J n jj* Aj (r) exp{z(>r • r)} df is orthogonal to %. 

The substitution of -V ? (A/| 5(f) \X. t ) + Vp A (A/| M(r) |Aj) for (A/| j(r) |Aj) in the 
magnetic scattering amplitude results in re- formulating £\ f ,Xi(^) in the eq. (|2.7p as 



-7n r 



2^ B ^' 



2 *A(A/ 



M(r) exp{i(i? • r)} (ir 



A, ) A x 



(2.11) 



M(r) is formally interpreted as the magnetic density operator since j\ \(f) — Vf A (A| M(r) |A) 
is understood as the stationary current density in the sample in the state |A). All the other 
contributions to the current density are necessarily irrotational, thus should be ascribed to 
external sources to comply with charge conservation. We can safely ignore them although 
they produce magnetic fields, for these might act solely on the neutron statistical ensemble. 
A magnetic field H.^(f) of dipolar nature is also created by the magnetic density. Within the 
sample it is equal to the opposite of the irrotational component of the magnetic density. 

(A| M(r) | A) behaves as the customary conceived magnetic dipole densities that emerge from 
the spatio-temporal coarse-graining of the Maxwell equations when the magnetic materials are 
approximated by continuous media. At the difference of these, by essence macroscopic but 
unequivocally defined, (A| M(r) |A) makes sense down to the subatomic scale but shows gauge 
invariance, since 



V r -A {(A| M(r) |A) + V P (A| *(f) |A)} = V P A (A| M(r) |A) 



(2.12) 



whatever the scalar field operator \I/(r) with constant quantum average at the boundaries of 
the sample and everywhere outside. This gauge freedom is a source of ambiguities for the 
reconstruction of the magnetic densities from the neutron data [3]. It is not eliminated by 
constraining the magnetic density to vanish outside the sample, where the solenoidal current 
density is null, nor by imposing that its volume integration gives the magnetic moment of the 
sample. It also survives over all the other electromagnetic equations for continuous media. 
As an example, the magnetic field created by (A| M(r) |A) + (A| ^(r) |A) is computed 
equal to H I ^ + ^^,(? ? ) = Hj^(r) — Vf (A| *f?(f) |A), which shows that the magnetic induction 
B(r) = fio{H(f) + (A| M(r) |A)} does not depend on the gauge choice for (A| M(r) |A). A 
gauge can be fixed with the equation \7? • (A| M(r) |A) = so as to retain only the solenoidal 
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component of the magnetic density, which is fully probed by the neutron, but then the concept 
is often not intuitive. Other gauges might reveal themselves more insightful depending on the 
physics under concern. 

Actually, the spin magnetism is naturally described by continuous distributions of point-like 
magnetic dipoles p,j at positions Fj 



(2.13) 



for these do exist objectively. All the other contributions to the magnetic density should be 
considered through the convection current density operator 



(2.14) 



of moving particles j of electric charge qj. Yj is the position operator and Vj the velocity 
operator of the j— particle. Vj = j?[H,Tj] = Vpv7i, where Ti. is the Hamiltonian of the sample. 

Nuclear magnetism is proportional to //n when electron magnetism is proportional to /zb, 
so can be disregarded. As a matter of fact, nuclear magnetism is rather probed through the 
neutron- nucleus interaction potential Vn- We shall henceforth consider only pure electron 
parameters : p,j = g c fiB Sj, where g c = 2.00232 w 2 is the electron gyromagnetic ratio, and 
qj = e Vj. Assuming further that the p— dependent terms in the Hamiltonian are dominated 
by the non relativistic kinetic energy, Vj = Accordingly, £x/,Ai(-^) m the eq. (|2.7[) ) writes 



1 



^[exp{i(if-?j)} 3j] 



Xi ) A x + 



+ 7nro m^^\ Xf 



[exp{i(>? • Tj)} p 3 + ft exp{i(x • Tj)}] 



(2.15) 



3 FORM FACTORS 



Let us consider the concrete instance where the sample is a collection of kinematically indepen- 
dent ions A p at the positions defined by the eigenvalues of R p . Any electron j then belongs to 
a single ion, which implies that any sample state |A^j) is the tensor product |A/,i) = (££) p (A^J 
of the ionic states |A^J and that (A/|Oj E a p |Aj) = (X^\Oj & A p \X^ ) whatever the one-electron 
operator 0,- e a associated to the j— electron of the ion A p . We deduce 



£\f,\i {%) = ~7o r o y^p * A -^ap,ap(x) A x\ exp{i(>?- R p )} 



(3.1) 



where 



M(r) exp{i(i< • r)} dr 



XI 



(3.2) 



defines the magnetic vector form factor of the ion A p . Using the spherical tensor operator 
formalism, it may be computed exactly, but this necessitates defining a number of mathematical 
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concepts and performing weighty algebraic manipulations, so will be detailed in the last section 
and in the appendix. In this section rather qualitative and simplified aspects will be discussed. 

Let T be the time inversion operator. T = Y K, where the unitary operator Y = 
exp{— i ir e y ■ J2j Sj} acts solely on the spin states and the antiunitary operator K of com- 
plex conjugation acts solely on the spatial states [1]. Let X be a pure imaginary operator, 
that is fully reversed by the time inversion: T X T _1 = —X. It follows that (-0/1 X = 
(,p f \ T _1 T X T^T = -<V>/| T- 1 X T \^) = -<T^/| X |T^)*, since, T being an- 
tiunitary, (ipf\T~ 1 x) — (V , /|1 ,+ x) = (Tipf\x)*- Assuming that X acts solely on the spatial 
states and that iV'/.il) are products of spatial states by spin states |C/,») or i if n °t, ex- 

panding over such basis states, the spin degrees of freedom can be ignored. We then may 
write (01 X |6> = -(K0| X |K&>* = -<K&| X |K&) = - f dr J ds &\r)* (r\ X \3)(3\Zf): 
If the Hamiltonian is Hcrmitian and if the states are non degenerate then the wave- 
functions (r|£/,») should be proportional to their complex conjugate, but K 2 = 1 so that 

= ex Pi il Pf,i}{r\^f,i)^ in which casc (01 X l&) = -exp{i(ip f - (Pi)}{£f\ X |&)*. If 
furthermore = £j then the matrix clement is an imaginary number. Using this, we may 
qualitatively infer that if the ionic states |A^) belong to a same electronic configuration (nl) 
then there should be no contribution to the matrix element of the convective current density 
from the radial degree of freedom of the electrons, since T jc(r) T _1 = — jc(^) but the di- 
agonal elements of jc(^) should be real measurable quantities. Explicit calculations confirm 



this by more precisely showing that if Vj 



then (nl\ jc{r) \nl) 



KM 



(l\ jc(rj |0, 



where R n i = (r\nl) is the radial wavefunction associated with the electronic configuration. If 
moreover the ionic states \X P t ) are in the same orbitally quenched state then for the same 
reasons there also should be no contribution to the matrix element of the convection cur- 
rent density from the angular degree of freedom of the electrons, in which case we may set 
M(f) = M s (r) = g c Mb V\ K^j - Sj so that 



E ex p{*(^- ? jp)} *j 



(3.3) 



Using the closure relation J2 V \v p )( r l p \ = 1 over the ionic states, this rewrites 



^exp{i(x-r jp )} J2\ V P )(V P \ *j 



jet 



Af = 



exp{i(i? • r 



Ji> 



)} \v P )(v v 



E 

jeAp 



(3.4) 



where the phase factor can be factorized because the matrix element of an one-electron operator 
over electrons states of a same configuration, which are linear combinations of antisymmetrized 
products of single electron orthogonal states, does not depend on which electron is chosen to 
compute it. A pure spatial state is insensitive to pure spin operators, therefore, because of 
the matrix element (?7 P | J2jeA p &j I^D' the closure relation |?7 P )(r] p \ = 1 gets limited to 
the ionic states rf with the same spatial component as the |A^). It then is inferred that 
{\ p \ exp{i(x- r jp )} |^>) ^ \\ p ) = \r,r). Accordingly 



(x) 



exp{i( 




i6A p 



AF 



(3.5) 
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where S p is the total spin operator of the ion A p . f\,{%) = (Ay | exp{j(i? • ?,•„)} |A^) defines 
the so called extracted scalar form factor. It is tempting to generalize this to non degenerate 
spatial states as 



(x) J M(f)exp{i(x. r)} dr A?\ = f p Xf ^) (^\ p f 



M(r) dr 



<> 



Af (3.6) 



which presupposes that the intra-atomic magnetic density M(r) is collinear. This often is a 
good approximation widely used in practice, in particular when the physical interest is rather 
focussed at the inter-atomic magnetic correlation, but should be considered with some caution. 
Non collinear magnetic densities are naturally expected in the presence of an antiferromagnctic 
order or any non collinear magnetic order, for there is no reason to believe that the magnetic 
orientation would abruptly change at well-defined inter-atomic boundaries. The non collinear 
character of the intra-atomic magnetic densities even come out in collinear ferromagnets or in 
paramagnets under magnetic field owing to the relativistic spin-orbit interaction which couples 
the magnetic densities to the local environments. Clear experimental evidences are revealed in 
precise measurements [4]. Strictly speaking, the concept of scalar form factors is limited. 

According to the exact computations detailed in the last section, the spherical components 
of the quantity {^j ^ ^ ^\ p f ,\ p {^) A >?} in the eq. (|3.1[) . for ionic states |Ay ■) belonging to the 
same electronic configuration (nl), can be put in the global form 

= £<?*(*)) £ c kK>M xP f \ p ) y<f (fi*) (3.7) 

' * K K>,Q< 

where Yq, are spherical harmonics over the spherical angles Ojj — (9%, ip%) of the scat- 

tering vector and C q K K , q<(A^, Af) gather matrix elements between angular momentum states 
and geometrical factors. At a first sight these appear algebraically complicated, involving coef- 
ficients of fractional parentage and summations over Wigner symbols, but in fact are formally 
not difficult to evaluate whatever the ionic states |Ajj). (jx {x)) are index K radial integrals : 

fe(x)) = / r 2 R 2 nl (r) 3K {nr) dr (R nl = (r\nl)) (3.8) 
Jo 

where R n i = (r\nl) is the radial wavefunction associated with the electronic configuration. 
These can be calculated at once and tabulated for each ion and valence by numerically solving 
the Hartree-Fock problem for the electrons interacting with the nucleus and between them 
in the one-electron central potential approximation. With the heavy ions the experimental 
measurements are enough precise to evidence significant deviation due to relativistic effects, 
for instance in the lanthanide or actinide series the f electrons are radially more expanded. 
Accordingly, the radial integrals for these ions are numerically calculated from the relativistic 
Dirac-Fock Hamiltonian [5,6]. This, in addition to the Dirac one-electron Hamiltonian and the 
electron-electron Coulomb interaction, includes the Breit interaction, which is treated as a first 
order perturbation, in order to partially take into account the relativistic character of the inter- 
action between the electrons : the Coulomb interaction is not Lorentz covariant but emerges as 
the leading term in an expansion of the interaction energy in powers of the fine structure con- 
stant obtained by the methods of quantum electrodynamics. As to be consistent the scattering 
amplitude also should be calculated on relativistic states from the relativistic expression of the 
current density operator j(r) = ■ c e atj 5(fj — f), where dj is the Dirac vector operator for 
the j— electron (notice that in the non relativistic limit this split back into the sum of a spin 
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and a convection current density operators) [1]. We shall not get into the details of the formu- 
lation, which is too weighty to be account for in these notes but can be found elsewhere [7,8]. 
An effective operator approach is considered in these works which allows expressing the matrix 
elements of relativistic operators over relativistic states as matrix elements of effective oper- 
ators over non relativistic states. It has been argued that weak corrections might occur in 
addition to the radial one, answerable primarily to the relativistic mass-correction term, which 
decreases the convection current density in the high-kinetic-energy region near the nucleus thus 
producing a more spatially extended current density, and secondarily to the relativistic spin- 
orbit term, which slightly changes the effective electronic g e factor. In practice, these effects 
are ignored and the relativistic scattering amplitude is obtained from the non relativistic one 
by merely replacing the radial integrals by their relativistic counterpart. 

When the inverse modulus x^ 1 of the scattering vector x is much larger than the mean 
radius of the radial wavefunction of the unpaired electrons giving rise to the magnetic scattering 
the magnetic vector form factor can be approximated by the very simplified expression 



(j (x))S + i[Oo(x)) + (i 2 (x))]L 



A' 



T(x) 



A 



(3.9) 



where S and L are the total spin and total orbital angular momentum operators of the ion. As 
discussed in the last section, this is inferred from the behavior of the spherical Bessel functions 
for small arguments and defines the dipole approximation to the scattering amplitude. Notice 
that the form factor then is isotropic in the reciprocal space, that is it depends solely on the 
modulus of the scattering vector and not on its spherical angles. If the unpaired electrons are in 
an orbital singlet state then T{>i) = (jo(x)) S, which gives a scalar form factor f(x) w (jo(x)} 
at small scattering angle. This often applies to the d electrons of the 3d transition metal ions, 
which, owing to their spatial extension, often experience strong crystalline electric field in the 
materials. Generally, the quenching is partially raised by the spin-orbit interaction. This acts 
as a perturbation weakly coupling the ground orbital state with excited orbital states, thus 
inducing a small orbital moment. The total magnetic moment of the ion then differs slightly 
from that of the spin contribution and comes out through a value of the electron gyromagnetic 
ratio g different from that of the naked electron g e . In order to keep the spin operator as 
the basic variable it is convenient so set (L + g e S) — gS, which determines g, in which case 
T{x) = i{(joM) g e S + [(jo(x)) + (j2(x))} L} is expressed in the form 



where 



T^x) = -gf(x) S 



/(x) = Oo(x)) + (l-| 



02M) 



(3.10) 



(3.11) 



Notice that we systematically set g e w 2. Unlike the d electrons of the 3d transition metal 
ions the f electrons of the lanthanide ions are much less extended, thus are subject to weaker 
crystalline electric field in the materials, whereas their spin-orbit interaction is strongly in- 
creased, becoming dominant. Accordingly, the total spin S and total orbital L moments are 
first coupled to form the angular moment J = S + L and the energy spectrum gets structured 
in terms of the multiplets \tLSJM) on which the crystalline electric field potential will act as 
a perturbation. Within each multiplct the Wigner-Eckart theorem allows writing 2S = gsJ, 
L = J and 2S + L = gjj, where g.j = gs + 5l, gs and gL are geometrically determined from 
g s P = 2S • J = S 2 + 2S • L = J 2 - L 2 + S 2 and g L P = LJ = L 2 + LS = Z+l^JL^ using 
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the identity J 2 = (S + L) 2 = S 2 + L 2 + 2S • L : 

g(g + l)-£(L + l) 1 £(L + 1)- 5(5+1) 

ffS = 1 + J(J+1) - 9i= 2 + 2J(J+1) (3 ' 12) 

Thus, if the initial |A?) = \rLSJMi) and final |A^) = \rLSJMf) states belong to the same 
multiplet (Jj = Jj = J) then 

^(*) = ^ffj/W J (3+3) 

where 

/(x) = <*>(*)> + ^0"2(*)> ( 3 - 14 ) 

If the initial |Af> = \rLS{JM)i) and final |A^> = \rLS{JM) f ) states belong to different 
multiplets (J; ^ J/) then (r J L5(JM) / |J|rL5( JM);) = so that (TLS(JM) f \S\TLS(JM)i) = 
-(TLS(JM) f \L\TLS(JM)i). It in this case follows that 

= \ f(JM) s .,(jM)A K ) S (3+5) 

where 

f(jM) f ,(jM) t (x) = CM*)) - CM*)) (3-16) 



So far the unpaired electrons giving rise to the magnetic scattering were assumed kine- 
matically independent from atomic center to atomic center, but this, in the real magnetic 
materials, is more the exception than the rule : the states of neighboring centers get mixed 
by the inter-atomic exchange interactions, there may be covalent transfer from ligand ions 
surrounding the magnetic ions, modifying the states of these and creating spin polarization on 
the ligand ions, the electrons under concern may tend to participate to metallic bonding, et 
cetera. The magnetic scattering then contains multi-center matrix elements and its compu- 
tation become more difficult. We shall not get into the details of the different instances, for 
these are quite varied, and shall only emphasize that the robustness of the one-center states 
is such that the use of the ionic-like form factors is often a good approximation. Thus, in 
most instances the intra-atomic exchange interactions are much larger than the inter-atomic 
exchange interactions, freezing out atomic or ionic angular momentum states on which these 
act as a perturbation, and the overlap mechanisms giving rise to these generally modify the 
spatial wavefunctions at large radial distance. Accordingly, the electrons preserve an atomic or 
ionic character and can be analyzed as if they were kinematically independent, at least at not 
too small modulus x of the scattering vector. Actually, even when more significant effects on 
the scattering amplitude might be expected, for instance in the presence of strong ion-ligand 
covalent transfer, the single ion analysis is worth performing for it might allow assessing the 
relevance of the multi-center contributions. The case of itinerant electrons must be treated 
separately, using Wannier functions. These show strong atomic-like character in the case of 
3d transition metals and 3d - 4f (or 5f) intermetallics and, ignoring the s electrons, spread 
essentially over the z nearest neighbors. Labeling these with vectors p, the magnetic scalar 
form factor fw{%) = /a(x){(1 — zW 2 + W 2 Y^p exp(? h ■ p)} can be assigned to each atomic 
center, where W is a Wannier expansion coefficient and /a(x) an ionic scalar form factor. 
The interesting point is that if i? is a reciprocal vector K then fw(K) = Ja(K), that is the 
itinerant nature will not be seen by elastic scattering [2]. In the case of 4d, 5d and 6d itinerant 
electrons, strong deviations from any ionic scalar form factor are generally found out even in 
the elastic scattering, signaling that the atomic integrity of these electrons then gets destroyed. 
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4 EXPERIMENTAL POLARIZATION RATIO METHOD 

Neutron experiments give access to the differential scattering cross section. This is given by 
the square modulus of the scattering amplitude A[{{kiVi)\i\ ~~> {(kfUf)X/)}], statistically 
averaged over the initial states of the sample and spin states of the incoming neutron and 
summed over all the final states of the sample and spin states of the outgoing neutron, 



— - — £ p^p»* z2 \ A 



dVLdEf h 

Ai,i/j \f,v f 



{(fci^)AJ •w {(kfff)Xf)} 



S 



f — +Ei- E f 

2m„ 1 



(4.1) 

where the Dirac delta generalized function expresses the constraint of energy conservation. The 
statistical weight p\ t of the initial states |A$) of the sample is generally given by the Boltzman 
factor 

_ expj-Ej/kBT} 

' J2 E ^M-Ei/k B T} [ ■ > 

The statistical weights p Ui of the spin states \vi) of the incoming neutron depends on the 
incoming neutron beam polarization. This is evidenced by building up the density operator 

P = £>iK>i| (4.3) 

Its matrix representative over basis vectors in the neutron spin state space is a 2 x 2 non 
singular matrix. Any such matrix is a linear combination of the unit 2x2 matrix and the 
Pauli matrices (cf. eq. (|C.16j) ). which altogether form a basis of the vector space of non singular 
2x2 matrices. Accordingly, the density operator may write p = itl + v • cr, where 1 is the unit 
operator in the neutron spin space and er the Pauli vector operator. The Cartesian components 
Cm (jn — x, y, z) of cf satisfy the commutation relations 

[ er m , cr„ J — 1 e cri, 

where e mnl — 1 if (to, n, I) is co-cyclic to (a;, y, z), = —1 if (m, n, I) is anti-cyclic to (x, y, z) and 
= otherwise, (cf. eq. (|C.9|) ) and the anticommutation relations 

Combining these we may also write <r m cr n = 6 m>n l + £ mnl(T i, which, since Tr[er m ] = 0, 
straightforwardly gives 

Tr[cr m cr n ] = 2S m , n and Tr[<x m <T„<x ; ] = 2i e mnl . 

Accordingly, u = Tr ^ and v = Tr ^ Tr[p<r], but Tr[p] — ^2 u .Pvt — 1 then u = \ and 
v = jTr[pa] 1 which precisely defines the polarization P of the incident neutron beam. Thus 

p=^l+Pa (4.4) 

The scattering amplitude, in view of the expressions of the different contributions to the 
neutron-sample interaction potential (cf. eq. (|2.2p , eq. (|2.3|) and eq. (|2.5p ). can always be 
put in the form 

A \{(ki^)Xi} ~» {(fc>/)A/)}l = {vf \a{Z)f ti 1 + 0(H) f>i ■ <t \ Vi ) (4.5) 
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where the quantities a(x)f t i and (3{>c)f^ are matrix elements of a scalar operator a(x) and a 
vector operator /3(x), which refer solely to the sample. Symbolizing a(x)l + /3(x) ■ <x by O, 
it follows that 

I>< E \ A - {(*>/)A/)}] f= E pMo + Wf)(vf\o\^) = 

= E< v /I° (E J = ?>[p C? + C7] (4.6) 

but, since Tr[a m ) = 0, Tr[cr m <x„] = 2£ m ,„ and Tr[<T m cr n cr / ] = 2i e" 1 ™', 

Tr[p O+O] = Tr[p {a(x)l + /3(x) • <r} + {c*(x)l + /3(x) • <?}] = a+(x)c*(x) + 

+ ■ + a+(x)(P • 0(H)) + (P ■ P + (x))a(Z) + » P • 03 + (x) A /3(x)) (4.7) 

by understanding that the trace is to be taken over the neutron spin states only. Using the 
identity 

Slhw + Ei- E f ] = — — / dt cxp{-i tut} exp{-i K f ' } (4.8) 

2ttTuj J_ 0o n 

(cf. eq. (|A.4p ). where huj = — ^ — ~> ^ hnally is inferred that 

-^^- = -L J dt exp{-*wt}|(a + (x)a(x,t)} + (/3 (x) ■ 0(H,t)J + 

+ (a + (x)(P • fe t))) + ((P • ^ + (x))a(x,i)) + i P . ((^ + (x) A £(*,<))) } (4.9) 

by defining a(x, i) = exp {^-jp} a(x) exp { } and p(x, i) = exp j 2 -^} p(x) exp { }> 
where 7i is the hamiltonian of the sample, and using the notation (• • • ) for the statistical aver- 
age over the initial states of the sample (summation over the final states is made by closure). 

The cross section for elastic scattering is obtained by taking the infinite time limit and that 
for coherent scattering by averaging over the nuclear isotope distributions and nuclear spin 
orientations. (a + (x)a:(x, oo)) = (o: + (x)) (a(x, oo)), since processes well separated in 
time get uncorrelated while (o:(x, oo)) = (a(x)) and (/3(x, oo)) = (/3(x)) under stationary 
conditions. If there is no net nuclear spin polarization, which amounts to assume that the 
temperature is not low enough to allow for nuclear spin order or that the hyperfine field is 
sufficiently weak so that the electronic magnetism does not polarize the nuclear spins or else 
if a magnetic field is present that this will not polarize the spins too importantly, then the 
nuclear scattering will contribute only to (a(x)) by the interaction term A p independent on 
the nuclear spins, averaged over the nuclear isotopes (cf. eq. 1|2.2|) ), In the case of a crystal 
with nuclei d at positions fd + R, where the vectors fd locate the nuclei in the unit cell and 
the vectors R — • n i^ L h ( n j e ^) define the lattice translations, we get 



(a(x)) = E ( E Ad ex P{* exp{i x-R}=Y J F N {K) 8{x - K) 



(4.10) 
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where K is a reciprocal lattice vector and Fn(K) defines the unit cell nuclear structure factor. 
The magnetic contribution to the elastic scattering, according to the eq. (|2.7p and the eq. (|2.1ip . 
is given by 



(pO*)) = ~7n r 



i(x ■ r)} df 



(4.11) 



In the case of a crystal with spatially periodic magnetic density, 



A> 



M(r) 



^^[/(^expiir-^^PA^Ai M(f-R) 



we get 



where 



(/3(x)) = ^y^U{r) 5{x-K + f) (x A F M (K - t) A x 



K 



(It, 5^ 



2^b 



PA; ( A 4 



M(r) exp{i(x • f)} df 



u.c. 



A, 



(4.12) 



(4.13) 



(4.14) 



defines the unit cell (U.C.) magnetic structure factor. If the magnetic periodicity is identical 
to the nuclear, U(f) = S(f), then the coherent elastic scattering cross section reads 



Fi;{K) 



Ft I (K) + F* N (K)(P-Fi; I (K)) + 



+ (P • Fh(K))*F N (K) + i P • (Ffr(K) A F^(K))} 5{x - K) (4.15) 



where F^{K) = K A F M (K) A K. By Fourier inverting F M {K), we find 



H^A, (A* 

Aj 



M(r) 



A , = 2 ^ B 1 J- F m (K) cxp{-t(K ■ f)} 
-7n ro Vu.c. - 

K 



(4.16) 



where Vu.c. is the unit cell volume, that is, up to gauge invariance, the magnetic density is 
fully determined from the knowledge of the scattering amplitude at the reciprocal lattice points 
only. 

The polarization dependent magnetic chiral contribution i P • {F^f (K) A F^ (K)) to the coher- 
ent elastic scattering cross section should be expected only for non collinear magnetic density 
without center of symmetry. If this case is excluded and if the F^(K) can be accurately 
determined from complementary measurements, for instance in the paramagnetic phase with- 
out magnetic field or by X-rays scattering, then the polarization dependent nuclear-magnetic 
interference contribution P • [F^(K)F^(K) + F N (K)F^j(K)*] can be used to allow for pre- 
cision measurements of the F^j(K) and thus of form factors. It indeed is an evidence that 
weak Fj^{K), amplified by Fn(K), will much more accurately be detected than the non chiral 
magnetic contribution F^*(K) ■ F^(K), especially at large scattering angle where the form 
factors tend to vanish out and where extreme sensitivity is often required to get insights about 
fine details of the magnetic density. A necessary condition of course is that the interference 
contribution itself does not cancel out, which is the case if Fn(K) and Fj^(K) are not in phase 
quadrature. The method applies to polarized paramagnets, to ferromagnets or ferrimagnets 
and to in-cell antiferromagnets with in-phase unit cell magnetic and nuclear structure factors. 
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Experimentally, the so-called polarization ratios R{K) between the scattered intensities for 
two opposite initial polarizations +P and — P of the incident neutron beam are measured on a 
series of reciprocal lattice points K (Bragg scattering) . The extraction from these data of the 
Fj-f{K) in general is not straightforward, for these are complex vector quantities, the direction, 
modulus and phase of which have to be determined, when only one real scalar quantity, R(K), 
is measured. Models then must be worked out providing with calculated Fm(K) and, from 
these, calculated R(K), to directly confront with the experimental values. A direct access to 
the Fm{K) nevertheless is possible in collinear magnets with a centric crystal structure. In 
this case, both Fn(K) and Fm(K) are real. If furthermore the polarization +P is set parallel 
to the magnetization vector then 

K = l + 2sin 2 a 7 (i?) + sinM 7 (i?))2 . = F N {K) 

l-2sin 2 a 7 (iv) + sin 2 a ( 7 (i?)) 2 ' F M {K) 

where a is the angle between the magnetization vector and the scattering vector K. This is 
easily solved for F M (K), which, it should be emphasized, is not the modulus of F M {K) but 
its algebraic value along the magnetization vector. In fact, we get two mathematical solutions, 
but, in practice, it most often is evident to recognize which of the two is physical and which 
is unphysical. The actual measurements arc performed under strong magnetic fields, which 
in the case of paramagnets allows inducing the largest magnetization and overcome eventual 
significant magnetic anisotropy and in the case of a sample with a net global magnetic moment 
allows avoiding depolarization by selecting a single magnetic domain and orient this along the 
magnetic field direction. The experimental corrections to be considered are essentially extinc- 
tion corrections, associated with amplitude coupling within perfect blocks (primary extinction) 
and between perfect blocks (secondary extinction), which, increasing proportionally to |fo;| -3 , 
is accurately estimated from measurements at different incident neutron wavelength 27r|fci| _1 . 
Corrections which are identical for the two incident neutron polarization (absorption, • • • ) are 
irrelevant and using filters the fractional wavelength contamination can be removed. Almost 
perfect incident neutron beam polarization can be produced whereas depolarization by the 
sample can be minimized by using samples elongated along the magnetization axis and cleanly 
polishing the surfaces. No description of specific experiments will be given in these notes, but 
overviews of investigations can be found in the literature [9] . Also of interest to consult is the 
methodology inspired from the spherical neutron polarimetry to measure the form factor in 
in-cell antiferromagnets with magnetic and nuclear structure factors in phase quadrature [10]. 

5 ALGEBRA OF IONIC FORM FACTORS 

An in-depth quantitative analysis is provided in this section of the magnetic vector form factor 
of atomic electrons, more precisely of the expression between the braces in the eq. (|3.ip . Calling 
to mind that x A [p, /(?)] = i?A {— ih Vp/(r)} is null for /(r) = exp{i(i?- r)}, this also writes 

S £A/ A p a p(x) Ax = 



E 



1 i 
exp{i(x ■ r jp )} ( — x A sj A X - pj 



A? (5.1) 



which distinguishes the contribution of the spin magnetic moments from the contribution of 
the convection current density, fj = Tj — R p is the position operator of an electron j of 
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A p with respect to the position of the nucleus of A p . Using the spherical tensor operator 
formalism [11,12], this quantity can be computed exactly [2,13-16]. The method basically 
is not complex but the associated algebra is unwieldy. A fortunate fact is that we have to 
deal only with matrix elements of sums of one-electron operators, which ultimately can be 
deduced from the matrix elements of these operators between uncoupled single-electron states 
= \{nlmism s y Si ). 

5.1 Single Electron Spin Vector form Factor 

Let us focus our attention on the spin vector form factor, which is associated with the operator 
X^jeA exp{i(i? • fj p )} Sj and examine the matrix element between uncoupled single-electron 
states 

& n r v r(£) = Unlmisrn s )f\ exp{i(x- r^)} Sj | {nlmism s )\\ (5.2) 

We shall temporarily drop the indexes p and j for notation convenience (notice that the index 
j is unnecessary for equivalent electrons since these cannot be distinguished from each other). 
The spatial wavefunction associated with \rj) — \nlmism s ) will be denoted (r \ nlmism s ) = 
(r IV | nlmism s ) = R n i(r) Y^(1V) \sm s ) = R n i(r) \lmism s ), where R n i(r) is the radial com- 
ponent and Y^Vlp) = (Hp | Imi) the angular component. Appropriately inserting the closure 
relation J \r)(r\ dr = J Q \r)(r\ dr J s ^ \np)(flp\ dflp = 1 and using the multipole expansion 

exp(i x ■ r) = 4TrJ2i K jK(*r) [1^ («*)] * Y£(n?) (5.3) 

KQ 

to separate the variables x = (x, fijj) = (x, 9%, tp%) and f = (r, Q?) = (r, 6?, tp?) (cf. eq. (|B.23[) ) 
the spherical components of ©„. .^(x) can be expressed in the form 

©^(x) = 4 7 r^i if 0>(x)) />i D#(n*)]* ((lm lS m s ) f | Yq s q | (Imism^) (5.4) 

KQ 

where 

/>oo 

(j*(*)>/,< = / r 2 R f nl (r) R l nl (r) j K {xr) dr (5.5) 



are dubbed index K radial integrals. The 2K + 1 operators Yq are the components of the 
spherical harmonic tensor operator Yk of order K. (Qp\ Yq = (Qp\ Yq (Tip) = 

^e P ,np Y Q (^r) by definition. s_i = -^{s x - i s y ), s = s z and s+i = -^(s x + « Sy) are 
the spherical components of the electron spin vector operator s. Applying the Wigner-Eckart 
theorem separately to Y# and to s, 

((lmism s ) f | Y£ s q | (l mi sm s ) t ) = (-) l f~ m 'f f K \ (l } \\ Y K \\ k) x 

V / (5.6) 

x{-)'f- m 't[ Sf Sl )(sf ||s || Si ) 

V —rn Sf q m Si J J 

(cf. Appendix ID.2[) with the help of the eq. (|D.13[) and the eq. (ID.15[) for the reduced matrix 
elements, we get 

6^, ni (x) = (e.l^^OxM)/,, [Y$(to a )]* (-r- m °f +mi f6 8f , Si [l,,K,h] x 

KQ (5 7) 

If K k\( l f K k \ ( s f 1 s t \ 



J V —TO/, Q mi i J \ —m s , q m 
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where [ll 1 , l e 2 2 , ■ ■ ■ ] is a standard abbreviation for {(2k + 1) £1 (2h + 1) £2 • • • } 2 • 6« / t?i (it) is null 
unless K satisfies the symmetry condition If + K + k = even integer on the 3jm symbol with 
all m null (cf. eq. (|C.57[) ) and the triangular condition K = If + U,lf + k — 1, ■ ■ ■ ,| If — h |. 

Owing to the behavior of the spherical Bessel functions for small arguments (cf. eq. (|B.24[) ) it 
is inferred that for it — > the expansion of 6^ Vi (H) over the spherical harmonics [Yq (£1%)]* 
should be dominated by the K = 0, Q = contribution, that is, 



[tfCW = (4.)- and ( _^ J ^ ) = (-)"-*, 8 lfM 8„ 



mi f ,mi. 



e n/.'?*( i? )^o ~ 0'oM)/,i <5 mi/ ,m i( (K)/ I s ? I (sm,)i) (5.8) 

which defines the dipole approximation to the spin vector form factor. Notice that this is null 
if If 7^ Zj. Generalization to multi-electron ions is straightforward. 

The spherical components &^ fr)l (ii) of the matrix element ^j?A 6^p ^p(x) A i? are given by 
©^(x) = & q Vf , m (Z) ~ ^^E(-) ? ^ S n^ i ( i? ) ( 5 - 9 ) 

where x g is the g— spherical component of the scattering vector. Calling to mind that Yq (CI&) = 
( cos 0s andY-^Qj?) = =p(J F )5(i)l sin 6^ expjiit^}, we may write n q = k (^-)iY^(Q z ) 
whence 



V 



x f ©-/ ni (x) = ( ^) F g 1 (0 if )^(-)«Fi(0 i? )6-% i (x) 



9 

K,Q 



9 

(4»r)*E( B )*H M ( J J o)(J J Q ( 5 - 10 ) 



(cf. eq. (|C.62j) ). The 3jm symbol with all m null is computed to — ^= for K = and to 

for X = 2 (cf. eq. I|C.76P ). It is null by symmetry for K = 1 (cf. eq. (|C5T|> ) and by triangular 

rr-n— — — 

conditions for any other K. It then follows, thanks to the eq. (|(J.68|) . since Y ° (fia) = (-L) 2 
and taking into account the zero-sum condition q + q + Q — on the second 3jm symbol, that 

i 

= few*) + h" 9 (f ) 2 ? ( J J - g - g ) 6 ^ ( * } (5 - n) 



Using for the 3jm symbol the algebraic form 



1 1 2 
q q -q-q 



(-)" 



9-9 



(2 + q + q)\(2-q-q)\ 



_30(l-g)!(l + 9 )!(l-g)!(l + g)!_ 
(cf. eq. (|C.67[) ) and for the spherical harmonics under concern the functional forms 



(5.12) 



i i 

' 15 \ 5 / 15 \ 5 

Y'l^Q,*) = T ( g- J sin 0* cos 0* exp{±i<p«}, F| 2 (fi») = f — j sin 2 5 exp{±2i</> 2 } 
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(cf. eq. (|B.7[) ), it finally is found that 
©^(x) = sin 2 6* + ^= sinfecosfe (e^6^ w (x) ~ e^S" 1 ^*)) (5.14) 

(*) = f 1 - 5 sin2 e ^(^) 

--Lsin^cos^ e^S^^x) + ^sin 2 S ^^©^(x) 
5.2 Single Electron Orbital Vector Form Factor 

No more difficulties arise when dealing with the orbital vector form factor associated with the 
operator £ ieAp exp{i(x • r ip )} (-^xApj) = £ i6A exp{i(x • r Jp )} (>(AV ? J] and, 

in particular, with the matrix element between uncoupled single-electron states 



(nlmisrris)? 



exp{i(x-?j )} 5AV f 



(nZmOO (5.15) 



but the computations are a little lengthier. A detailed description of these when the states 
r]f and rji belong to a same electronic configuration, that is for (nl)f = (nV)i, is provided at 
several places in the literature [2,14-16] with only faint differences in the formulations. Let us 
replicate one of the approaches and re-write the matrix element 



nlmu 



1 



exp{i(x • r)} (x A V?) 



nlmi 



(5.16) 



- — j ^(nlmij, exp{i(x-r)} (xAV?) nlmi^j — (nlm^ (xAV?) exp{— i(x ■ f)} nlmi^j j 



1 



dr exp{i(x- ?)} i? 2 ; (r) x 



x M* ((* A V?) Yl h (IV)) - 5^ (IV) ((x A V?) (<V)]*) } 

The indexes p and j are temporarily dropped and the spin quantum numbers temporarily 
ignored for notation convenience. The second equality gets obvious when calling to mind 
that Vr) f VO {ijf \0\ra) = {Vi\0 + \Vf)* and VO x V0 2 (OiO a )+ = 2 + Oi + , where 0+ 
symbolizes the adjoint of O and that [exp{i(x • ?)} x A V?] + = — x A V? exp{— i(x ■ f)}, 
since (exp{i(x • r)}) + = exp{— i[x ■ r} and (V?) + = — V?. The third equality is derived by 
appropriate insertion of the closure relation J*|r)(r| df = J Q \r)(r\ dr J s ^ |IV)(IV| dQ? = 1. 

Using the spherical components of the vector operator x A Vp, 



(x A v ? ) 9 - Y, 



91,92 



1 1 1 

9i <?2 -Q 



(cf. eq. (|D.7[> ) and the multipole expansion 



exp(i x-r) = 4mY,i K 3K{xr) [y<f (n a )] * 1^ (IV) 

KQ 



(5.17) 



(5.18) 
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to separate the variables i? = (ft, 0,%) — (x, 6%, tp%) and r = (r, fV) = (r, 6^, ^) (cf. cq. (|B.23[) ). 
the spherical components of £„^ m( (i?) can be expressed in the form 



„„ i «(-) 9 \/6 



111 

91 92 -9 



(5.19) 



KQ 



x y dfiK(xr) i£,(r) ^(^{[i^/fV)]* (v 92 Y l mi xn P ))-Y^\n P ) (v 92 ^(IV)]*)} 

As from the identity 



V =-\Vl r 



V? = - d r (d r r) + 4 V 2 



92 2 
2 "^L - 9 r (9 r r) 



(5.20) 



I sin &■* v 



r A \ sin ( 



and since spherical harmonics are eigenstates of the angular laplacian V 2 , „ (cf. Appendix IB")) 



vg lV , ^(rv) = Ki + i)^(JV) (5.2i) 

we may write 

(K^OV)]* (v 92 Y l mi xn F ))-Y^in P ) (v 92 Pl./JV)]*)} = 

5 { K lf (Of)]* (V| [r 92 y^ (IV)] ) - Y^ h (IV) (v 2 - [r q2 [Y^ (IV)]*] ) } (5.22) 

but, substituting for r ?2 its expression r (j§-) 2 Y q2 (flp) in spherical coordinates (r, <pf) and 
taking into account the identities 



-in 



i i ? - i 



92 m q 



i / i + i 

q 2 m q 



(5.23) 



(cf. eq. (fCT62l) and eq. (jCJ6l) ) and 



V 2 - (r y_\(IV)) 



it is inferred that 



1 a , a , fc(fc + l) 



(«V)]* (v| [r g2 y^ (iv)] ) - y^ (IV) (v 2 ~ 



) (r Y_yiV)) = ~(fc - l)(fc + 2) y* 9 (IV) (5.24) 

[YL lf (n P )r])} = k-) l+1+m 'f x 



l*(i-2)(l + l) Y!_-'(Q f ) 

- {i + i)k(i + 3) yit^iv) 



i i 1-1 

92 mii q 



Y'm 



iHi-2)(i + i) y^(iv) 

-(Z + l)% + 3) Y 1 ^ 1 ^) 



1 I l + l 

92 m h q 

1 I l-l 

92 -mif 9 

1 / z + i 

92 g 



(5.25) 
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The angular part of the integral over r in £L, (it) thus is reduced to integrals of products 
of three spherical harmonics (cf. eq. (|C.63p ) and therefore effortlessly found out to be 

<m f y^(o r ) {[Yi if (a?)]* (v 92 (fi r -)) - (n,) (v 92 OV)]*) } = 



x J [Z(2Z-1)]'(Z-2)(Z + 1) 



X J Z-1 




if Z 1 Z i-l\_/JT Z Z-l\/l I Z-1 

Q -mif -q J \ ?2 mii 1 J \ Q ™/, ~3 / \ ?2 -mj, Q 



[(Z + l)(2Z + 3)] 3 Z(Z + 3) 



if Z Z+l 




if Z Z + 1 \ / 1 Z Z + 1 \ _ / K I Z + lWl Z Z + I 
Q -mi f -q J \ qi q J \ Q raj, -q J \ q 2 -mi f q 

(5.26) 

K must be an odd integer, owing to the the symmetry condition K + Z + Z — 1 = even integer 
on the 3jm symbol with all m null (cf. eq. (|C.57[) ). Using the identity 

\p,_Vq( I 1 l±l \ ( K I l±l 

4^ ' \ x q 2 q ) \ Q y —q 

(cf. eq. (IC.81j) ). where the triangular condition for the non vanishing of the 3jm symbols 
imposes that K ' = K , K ± 1, then finding out that K' must be an odd integer because of the 
factorization 

Z Z K' \ ( I I K' \ ,„ , ,2i+K's ( I IK' 



n> I \ < )' I -(±-(-r +K ) A, ) (5-28) 

so that K' — K only, and finally taking advantage of the algebraic expression for the 3jm sym- 
bols with all m null (cf. eq. (|C.76j) ) and the following formula for the 6j symbol (cf. eq. (|C.89|) ) 



Z K z 
1 z-1 K 



l+K+z 



2(1 + K + z + 1)(K + z - l)(l + z - K)(l + K-Z+1) 



2K(2K + 1)(2K + 2)(2z - l)2z(2z + 1) 

(5.29) 

to notice that 



2t(if,Z) = [1(21- l)] 5 



K I Z-1 \ J Z K I 
1 Z-1 K 



= [(/ + 1M2/ + M.P ( ^ o ' o 1 ) { 1 f ' <■••*» 
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we get the more compact expression 

dn P y* (jv) {\Ymi f (*V)]* (v ga Yi u (iv)) - (IV) (v ga [Y^ (IV)]*) } = 

2 / 1 V [^M 8 ] R W 2W)£(-) ' (Wm, I luii £ ) ( g * ) (5.31) 



r V 47r 



The radial part of the integration over r in the expression of £^r~ m (i?) is more straightfor- 
wardly found out to be 

dr -Jk^t) R 2 nl (r) = —p— {<j>-i(*)) + (5.32) 
r 2K + 1 



where 

/>00 

O'atW) = / * (5-33) 



thanks to a recursion relation for the spherical Bessel functions (cf. eq. (|B.25|) V Gathering the 
equations and using the multipole expansion 

£ * ^ (5-34) 
(cf. eq. (|C.62p ) and the identity 

v ( K 1 M f 1 * ^V 1 1 1 

' \ Q 92 -Q' J\ qi -Q -Q J \ Ql qi -q 

, K f K K 1 





A' 




(J 








(cf. eq. (|C.82|| ) it finally is inferred, exploiting the symmetries and triangle conditions of the 
involved 3jm symbols, that 

£^, mi< (x) = (8n)H-) l+1 [l 2 ] El E ^ +1 {<i*-iM) + fei(x))} [^ 2 ] B(#,0 x 

K.Q I *,Q' 





K 






if 




(i 





o) 


n 


1 


i} 



(KQ'lmn | Zm;/) (KQKQ' \ lq) j [if] lj(n») (5.36) 

if = if ± 1 must be an even integer. 
When vc. — > the expansion may be limited to the lowest K = 1 order, in which case K = 0, 2 
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and 



(i?)^ Q = (8^(-) i+1 W] i*{{j {x)) + (j 2 (x)» 3 11(1,1) x 



1 1 




{111 }l>QW I W> (ooiQ' I ig> 



J o o ) { i i i f E (1Q ' Zmii 1 ' m < /} ^ 2Q1Q ' 1 ^ ^ Y i^ 



(5.37) 



It is a matter of elementary algebra to compute 



st(i,0 = VK2f-i) 



1 o 1 




(001Q' | lg) = 



1 

7! 



1 


i 


I 














o 1 ) 


{ 


1 





1 




1 


1 


1 


1 


}- 


~3 


r 


U) 






-Q' 



I 1 I 
1 l-l 1 

1 2 1 




(-) 



'+i [7-21 



6 

1 2 1 
1 1 1 



(2 + Q'- 9 )!(2-Q' + g)! 



10(1 -Q')!(l + Q')!(l -«)!(! + 9)! 



1 

" 6 

Iq-Q' 



taking advantage of the symmetry properties given in the eqs. (|C.55P - (|C.56[) and (|C83|) - (|C84|) 
and using the eqs. (|C.67jl . (|C.76jl and (|C.89|1 to get the algebraic expression for the required 3jm 
symbols and 6j symbols as well as for the quantity 21(1, 1). Considering then the functional form 
of the spherical harmonics displayed in the eq. (|5.13| , and calling to mind the Wigner-Eckart 
theorem to show, with the help of the eq. (|D.14[) for the reduced matrix element, that 



(lQ'lmu | Imif) 



(lmif \ Iq> \lmu 



(5.38) 



where 1 is the spatial angular momentum operator of the electron, we get, by comparison with 
the eq. 1(5TT4)) . 



{O'oO)) + (.h(x))}( lmif \ -^xMAm: \lmu ) 



(5.39) 



which defines the dipole approximation to the orbital vector form factor. The generalization 
of this approximation to more than one electron is straightforward and gathering this with 
the dipole approximation for the spin vector form factor immediately leads to the eq. (|3.9p . 
This often is used for an estimation of the orbital contribution to the experimentally measured 
magnetic density in magnetic materials and is the origin of the so-called (jo) — (72) analysis. 
An exact computation, for instance for the lanthanide ions in the ground multiplet level of 
maximum azimuthal quantum number, shows that the obtained values may roughly be good 
in some instances, as for Nd 3+ or Pr 3+ , but not always, as for Dy 3+ or Ho 3+ , owing to the 
additional contributions to the coefficient of (J2), from the quadrupolar term (K = 2) in the 
expansion of &^ ^(x) and the octupolar term (K = 3, K = 2) in the expansion of £ j; ^ j;P (x). 



5.3 Multi Electron Magnetic Vector form Factor 

Let now us generalize to correlated electrons within a same electronic configuration, that is to 
equivalent electrons. Adopting the = \(vLMlSMs)^ ,) quantization scheme, associated 
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with the total orbital L = J^j lj an d total spin S = Ylj &j kinetic moments of the equivalent 
electrons, we are led to compute the matrix elements 



x A ^p, s p(>?) A x = UvLMlSMs)) {vLM l SM s )\ 



p\ - 



(vLM L SM s ) p f 



E 

i6A p 



exp{i(£ -rj )} 



~ X ASj A X — r-rX A p 7 

2 fr>< 2 



t> groups additional quantum numbers, to be precise irreducible representations of certain 
groups, discriminating between the states with identical LM^SMs quantum numbers. The 
operator F 1 - is the sum F 1 - = Y^jLi °f one-electron operators. Using the concept of 
fractional parentage coefficients that emerge when building up a correlated n-electron state 
from correlated (n-l)-electron states and single electron states, its matrix elements may be 
calculated from the single electron matrix elements. The methodology details, too lenghty are 
provided in the Appendix [D] and the Appendix [EJ Use must be made of the eq. (|E.6|) for 
the convection current density contribution since we have to consider a tensor operator acting 
only on the orbital state space and of the eq. (1E.9|) for the spin current density contribution 
because then we have to deal with double tensors, which behaves as a rank 1 tensor in the 
spin space and as a rank K tensor in the orbital space owing to the spherical harmonic tensor 
that emerge from the mulipole expansion of the exponential operator to separate the spatial 
electron variables from the scattering vector parameter. The spin contribution reads 

{{vLM L SM s f f \ ]T exp{ l (x.f, p )} ( Sj ) q \(vLM L SM s )») = (6rr)* [Z 2 ] x 

jeA p 



KQ 



I K I 





■Af£(0/{|0)(0|}0i)(-) S+L [SfSiLjU 



Sf 1 Si 

-M Sf g M Si 

Sf 1 Si 

1 c i 

2 u 2 



K L; 



-M Lf Q M u 



U K L f 
I L f I 



(5.40) 



while the orbital contribution 



[vLM L SM s f f 



E 



3 et 



hx 2 



K A pj ) 



{vLMlSMs)? 



= [i 3 L f u] s Sf ,sMsf,M S A-) l - Lf+Li - MLi x 

E ( E lK+1 (OV-iW) + [K 2 ] *{KM-r K+Q ' (KQKQ' I Iq 

K.Q I K,Q> 





K 






A" 




(1 





1) 


K 


1 


i} 



Li A L/ 
Mm Q' -M L/ 

Li K L f \\ 
I L I | | 



A 



y|(^) (5.4i) 



These matrix elements can be computed in the \(vLSJMj)^ ^) quantization scheme, either 
directly or through the base transformation 



Mi M s 



(5.42) 
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What matters at this stage is actually not the exact and full expression of the different matrix 
elements in the different quantization schemes, as these can be found at several places in the 
literature [2,14-16], but rather to get aware that these essentially depend on the single electron 
matrix elements and the few coefficients of fractional parentage. Although the calculation of 
these coefficients may not be easy, it has, at any rate, to be done only once. Finally, it is in 
principle not difficult to generalize the computations to the case where the scattering event 
induces inter-configuration electron transfer. 

APPENDIX 



A FOURIER TRANSFORMS 

Let T( fc ) be a tensor field of rank k over the 3-D real space R 3 with integrable components 
T[ fc ^ (I = —k, ~k + 1, ■ ■ ■ , k). Its Fourier transform T{T^} is the tensor field of rank k over 
R 3 — in fact over the dual of R 3 but identified with R 3 — the components ^{T^ji (I = 
—k, — k + 1, • • • , k) of which are the Fourier transforms .F{T^} of the , that is defined by 

T{T\ k) }(q)= f T ; (fe) (r) cxp{-t(q-r)} dr (A.l) 
Jm 3 

TPV) = 77^3/ exp{itf - r)} dq 

(A) J r3 

is angular independent in the dual space, F{T9^}(q) = J 7 {T^}(q = \q\) \/q, if and 
only if Tj fc ^ is angular independent in the direct space, T[ k \r) — T[ k \r = | r |) Vr, in which 
case 

A-TT /*°° 1 /*°° 

;F{T< fc >}(«) = — / rT< fc >(r)sin(«r) dr | T ; (fc) (r) - _ / g^{lf >}(«) sinfar) dq. 
1 Jo Zir r J 

As an example, 

JF{i}(<7) = lim JF{^-^ — — —}(q) = — lim f exp(— er) sm(qr) dr = 
r e^o r q e^o J Q 

4-7T f°° exp(iqr) — cxp(-iqr) An 1 1 1 An 

= — hm / exp(— er) dr = — lim ( — = — , 

q e^o J 2i q e^o e - iq e + iq 2i q 2 



that is 



n\m = ^ 1 l^l^^m-mdq ( a. 2) 



S, • • • symbolizing scalar fields (k = 0), it is an easy matter to demonstrate that 

1 

(2*) 

where 



* S} = F{$}F{E} T{1>~} = 77J -^jr{$} * jr{ S } (A.3) 



($*£)(r)= / $(r-s)S(s) ds | * = f F{<f>}(q - k)F{E}(k) dd 
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The neutral element for the field convolution * is the Dirac generalized function in the R 3 space 
since by definition \f\ Jr S{r — s)x(s) ds = that is V% S * x — X- We deduce 

T{5{r)} = \ and = * ^{1} = _L_ /" exp{-i(g • r)} dr. (A.4) 

(2tt) j {2ttY J r3 



G, H, • • • symbolizing vector fields (k — 1) with components G/, H/, •••(/ = —1, 0, 1), 



:r\ <i> g } = :r\ <i> J ;r{ g j ;/-{ <t»g ) = tA^W * ^{G} 



^{^ G ; * H;} = ^{G} • ^{H} T{G • H} = -i- ^^{ Gi } * ^{H ; } 

1 \ f 1 



(A.5) 



J2 t lmn HGi * H m }e„ = ^{G} A ^{H} ^{G A H} = ^3 J2 e lmn HGi} * ^{H m }e 



where e' mn = 1 if (I, m, n) is co-cyclic to (1, 2,3), = —1 if (I, m, n) is anti-cyclic to (1, 2, 3) and 
= otherwise. 

9 rg Tp^(r) symbolizing the derivation of T^ k \r) with respect to rg = e • r (e.g. = x,y,z) 
and assuming that Tf (r) tends algebraically to zero at infinity, 



HdrM k) }(q) 



d r M k \r)] exp{-i(<f ■ f)} df 



T^ k \f) exp{— i(q ■ r)} 



T ; (fc) (f) [d r , exp{-i(q- f)}} dr = (ie-q) T{T\ K, }{q) 



(fe)i 



We find in particular that 

$} = i qT{1>}, F{A? $} = -q 2 T{<£}, 
■ G} = i q ■ T{G}, A G} = i q A ^{G} 



(A.6) 



These formulas are helpful for solving a number of partial differential equations. An example 
is A? = — 47r5(r), which is immediately solved as $(r) = -. Another example is the pair 
of maxwell equations 

V? A B(r) = fMoJ(f) V P • B(r) = (A.7) 

which must be solved to find the magnetic induction B(r) created by a time independent 
current density j (r) . We get in the dual space 

i q A T{B} = Mo ^{j) i g ■ ^{B} = 0. 

(cf. eq. (|A.6jl ) The second equation, which merely is gauge fixing for B(r), tells that there always 
exists a vector field A(r) such that .F{B} = i gA J-{A}, or equivalently B(r) = V-a A(r). The 
first equation then gives 

i q A (i <f A ^{A}) = (i q- T{A}) 1 q + q 2 T{A} = Mo^{j), 
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since V(w, v, w) u A (v A w) = (u ■ w)v — (u ■ v)w. A(f) by definition is subject to gauge 
invariance and we are free to fix this as i q ■ ^{A} = 0, or equivalently Vf • A(r) = 0. We 
deduce 

q z 47T r 

H B} = £° i gAf{V{j) = ^ * Q *j> - ^ ^{V F A [ 1*1]} 

47r r 47T r 47T r 

(cf. eq. (|A.2|) . eq. (|A.5|) and eq. (|A.6|) ). that is 



/,!! ; m^^^-ds (a. 



B (rl = ^v,A / 

47T 

Any vector u is the unique sum of its longitudinal component -^(q-u) q and its transverse 
component ^ (<f A tt A q) over the vector q. On transposing this to the Fourier transform of a 
vector field G, we find 

F{G} = \[(q- T{G}) q + (<f A ^{G} A g)] 

= i- 9 (i 9 ' -^{G}) + (i q A (i <f A ^{G}))] 

47r r 

= -i <f (A * V s - • G}) + i q A (-L ^{i * A G}) 

47r r 47T r 

= (^- - * Vy • G) + V r - A (-j- - * Vj- A G)}, 

47T r 4-7T r 

that is G is the unique sum G = Gj + Gs of an irrotational component (Vp A G/(r) = 0) and 
a solenoidal component (Vf • Gs(r) = 0) 

G(r) - G/(f) + G s (r) = -V F T(r) + V? A S(r) (A.9) 

where 

4?r J R 3 I r - s I 4?r 7 R 3 | r - s | 

T is determined up to a real constant and S up to the gradient of an arbitrary scalar field. 

B MULTIPOLE EXPANSIONS 

A multipole expansion by definition is performed over the collection {3^ i }_« TO <; ) 1=0,1.2, ■■■ of 
spherical harmonics and by essence is substantiated from the Sturm-Liouville spectral problem 

{V^+A}/(0,y>) = O (B.l) 

on the unit 2-sphere £2 with the boundary condition that / must be finite everywhere on S2. 
This belongs to a wide class of mathematically solved problems, but can also be approached 
more intuitively by interpreting the operator 

-V^ = - (-^de (sm9 d e ) + -±-d% 
\smff smO z * 
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as the square L of an orbital moment L and the operators 

f cos 6 1 

L z = -idtp and L± = ±exp(±i(p) <d g ± i-r-gdp > (B.2) 

as the components of L on the polar axis (z) of S2 and on the two associated orthogonal helicity 
axes (±). [V^,LJ = 0, which implies that 3 Y£ : V^Y* = _ Ay A and ^yA = m yA 
[L z , L±] = ±L ± so that L z (L ± Y r ^) = (m±l)(L ± Y^), which suggests that L± are indeed ladder 
operators. [V 2 , , L±] = 0, which tells that Y£ and h±Y£ belong to the same eigenspace labeled 

by A. [L+,L_j = 2L Z and -V 2 ,^ = ±(L+L_ + L_L+) + L 2 then L T L± = -V 2 .^ - L 2 T L 2 
and 

< f sintf d<? f *d^ [L ± ^(0,^)]*L ± r^(0^) = (L±F^|L ± y^) - (y^l^y*) = 
Jo Jo 

= <l^|-V^-L^L,|y^) = (A-m(m±l)) / sin0 d0 / d^ [Y*(0, p)]* Y^(0, p), 

Jo Jo 

which indicates that if is normalized then so are Y^ ±1 and that once A is fixed m is 
bounded. If I is the upper positive bound then (L+Y^JL+Y^^) = so A = 1(1 + 1) and the 
lower negative bound is necessarily —I, hence the more appropriate notation Y^ (—1 < m < I). 
An outcome is that V 3 p £ N 3 ? £ N : ^+ Y m = and L q _ Y^ = 0, that is m + p = I and 
m — q — —I, which implies that 2m = p- geN and 21 = p + q G N. Since L z Y^ = —id v Y r l n — 
mY^ the functional form of Y^ should be Y^(0, <p) — A l m (9) exp(imy) . Accordingly, if we 
further require that Y^ must be a 2tt— periodic function in the angle <p then m itself must be 
an integer and so 1 as well : the allowed values of A are 1(1 + 1), I = 0, 1,2,- • • and for each I the 
allowed values of m are — I < m = 0, ±1, ±2, ■•■</. Negative 1 are irrelevant because Z(Z + 1) is 
invariant under the I — ► — (Z + 1) transformation. Each (7, ro) distinguishes a single Y^, which 
means that the operators V# and L z form a complete set of commuting operators. Now from 
L±Y|, = it is inferred that d e A l ±l (9) = l^A l ±l (0) so Y^ =±l (e,tp) = C ±i sin 1 6 exp(±il<p), 
where \C±i | 2 is fixed from the normalization condition (Y_y | Y± t ) = 27r|C±; | 2 sin 2/+1 8 d6 = 1 
and C±; from the phase convention C; = (— ) |C±/|. It is then advantageous to make use of 
the identity 

2 f 2 cos 2zi+1 9 sin 2z2+1 9 d6= - — Zl! ^ — V z 7 - G C : ftfo) > -1 (j = 1, 2) (B.3) 
7 (zi + 22 + 1)! J J 

which is easily deduced from the double integral x^y^ exp(—(x 2 + y 2 )) dxdy by per- 

forming at first the (u = x 2 , v = y 2 ) variable change then the (x — rcos#, y = rshi9) 
variable change : we are led to the equality | J °° u^ - exp(— it) dit | / °° d^" exp(— u) = 
Jo°° exp(— r 2 ) dr J Q 2 cos^ #sin^ 6 dO from which the identity immediately follows, with 

£ = 2z\ + 1 and £ = 2z2 + 1, by calling to mind that z 1— > Jp 00 i z exp(— £) di = z! defines 
the factorial function : t z = ^(^ cx P 9( z ) io g (t) and |J o °°^exp(-t) dt| < J °° exp(-i) d* 
so J °° rj z exp(— rj) dt makes sense for all complex numbers z with real part JJ(z) > — 1 while 
integrating by parts the functional relation z! — z(z — 1)! is recovered and in fact might be 
used to extend by analytic continuation the concept of factorial to the whole field C of complex 
numbers. When z\ = —Z2 = z with —1 < 5R(z) < 1, the identity (IB.3P reads 

z\(-z)\= 2 (cot 2 BY d(sm 2 6) = / -dw=l 

Jo Jo (l+w) 2 
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where w = cot 2 6. 1 can be evaluated by integration over the complex plane : the analytic 
function s e C i— > g(s) = admits —1 as a pole of second order and as a branching 

point. g(s) is made uniform by cutting the complex plane along the positive real axis so that 
< arg(s) < 2tt and s z = \s\ z exp{iz arg(s)} with arg(s) = + at the upper adherence to the 
cut. Let £> be the domain formed by a disk of radius R centered about s — from which the 
cut along the positive real axis and a small disk of radius r centered about s = — 1 are removed 
and let £ be the external and c the internal closed contours limiting 2), oriented anticlockwise, 
then 

Vi? V r < R J g{s) ds = J g{s) ds 
but / g(s) ds - R ~ > °°> X (1 — exp{i 2tt z}) and — — / g(s) ds r-> °> (exp{i ir z})(— z) 

J<t l 27T J c 

since \g(s)\ « {s^- 2 < ^ (|9ft(z)| < 1) when ^ oo and 5 (.s = -1+e) = exp{i tt z}^-^ = 
exp{i 7r z}(^ — | + ■ • • ) when r — > 0. As a result 

z | = (B.4) 

sm ttz 

true for all z in C by analytic continuation. We in particular have f |j ! (— i) ! = 5, which implies 

that (-£)! = 0F and (£)! = ^ since (i)! = ± (-§)! and (-±)! > 0. When Zl = z 2 = z 
with 5R(z) > — 1, the identity (IB.3|) reads 



2 / (sin6>cos6>r z+i d<9 = ' = 2-2— r / (sin<^)^ +i 



o (2z + l)! ^+^0 v r; r 2 2z + 1 (z + i)! 

from which it is inferred that 

V¥(2z + 1)! = 2 2z+1 z\ ^ + ^ ! or V^(2^)! 2 2z z\ (z - 0! 

true for all z in C by analytic continuation. 



(B.5) 



Let us focus our attention back to the functional form of the spherical harmonics. Up to phase 
convention 

Li^m = m + 1) - Mm ± l)]*YUi = [(I T m)(l ±m + 1)]*F^ ±1 
since L z (L ± Yi) = (m± 1)(L ± ^) and (L^L^) = (f(Z + 1) - m(m±l){Yj n \Yj n ), while 

L±1^(0,V?) =±exp{±i^}(9 e ±i^^^) A l m {9)exp[im<p) = 

= ±exp{±^}(<9 e =Fm^r) A^(0) exp(irrup) = ± exp{±^}(sin ±m «9 e sin^™ 0) 1^(0, y>) 
sm p 

Accordingly, 

p times 



'((±m)!((^m-j))! 
(Z =F m)\(l ±m+p)\ 



{(Tysm^vhd^sm^eeM+iP^Yl&ip) (B.6) 
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where d cos6 = j^d e = -^d e . We deduce Y*{0, = -^^Jm^ d[ os9 ^ 1 9 and 
Yp(6, ip) = \J (j^pji (—) p sm p 9 d p osg exjp{ipp} Yq(0,<p). A more standard expression is 



121+1(1-771)., 



where 



YM0,<p) = (-) m J ^ ^ ^ i^( C0S g)exp(imy) (B.7) 



are dubbed order m associated functions of the Legendre polynomial Pg of degree I. 



I 



-/ + m)!(/!) 2 ( a ;-l)'- fc (x + l)'- fc +'" (fe ^+ m) 2 _ (/ + m)! 
jfe!(i + m-fc)!(i-A!)!(fe-m)! 1 ' (I — m)\ x 



,_ () jfe!(i + m-fc)!(i-A!)!(fe-m)! ~ k j (i-m)! x 1 

thus 

= (-rfr^f^^) i [^ l (c^)]* = (-r^ m (c,V') (b.q) 

Argument symmetry properties : Y^(C,-V>) = [Y^ (C, V)] * , ^m(-CV') = (-) m ^m(C», 
^(-C,-V>) - Yl m (C,i/>), Yi(Tr-C^) = (-) l+m Yi(C^), ^(Ctt + V) - H m ll(C». 

Vg and L z are self-ajoint operators, which implies that the Y^ are orthonormal, namely 

/ sm9d9 dp [Y l l (9,p)]*Y q k (9,<p)=5 lk 6 mq (B.10) 
Jo Jo 

Let £| 2 be the vector space of square integrable functions on S2, £| (N) the subspace engen- 
dered by (=0,1,2,... ,n and ^(N) 1 - its supplement : £| 2 ^ £| 2 (N) © C% 2 (N)^. 

N / 

V ff e £| 2 = .g - E E < r ™ Iff) e £| 2 (N)- 1 - 

(=0 m=-i 

so ^l-v^ M s . nce m . n (/l-Jp) =(N + 1)(N + 2) 

(/in|/in) /e£ 2 S2 (N)^ (/|/) 

N I 

but (fc N | - V 2 >ip h N ) = (g\ - Vj^g) - E E \( Y i\9)\ 2 Kl + 1) < (<?l - V^ ff > 

;=0 ro=-i 

then, (5I - V e 5) being N-independcnt, < lim (/i N |/i N ) < lim — A y — = 0, 

N-+00 N^oo y IN + 1 ){ IN + Z) 

which establishes the completeness of the semi-infinite set {Y^}_/< m <^ ;=o,i,2,-.. hi £| to be 
precise 

V 5 e£| 2 <7(0,^) = EE / sinCdC/ # [ll(C»]* 5 (C» hK^V) (B.11) 
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An equivalent form of the completeness relation which avoids specifying any function g in C s ^ 



is 



E K(C»r Yl(e,<p) = -^~5(9-0~5(¥>-i>) (B-12) 
1=0 m=-l 

where S is the Dirac generalized function defined on the < 9, C < 71 an d < ip, ip < 2tt 
compact domains. Notice that -^-gS(9 — Q = 6(cos9 — cosC) and more generally if / is a 
uniform function with derivative d x f and if {cii} are the roots of the equation f(x) = then 
= |a j(a )\ ( x ~ a *)' m com pact domains T>\ of the real space R and i5 in compact 
domains D3 of the 3-D real space R 3 are defined such that V/ J v S(x — u)f(x) dx — f(u) and 
V$ J v 5(r — s)$(r) dr = If (r,9,ip) are the spherical coordinates of r and (s, £,-)/>) those 

of s then dr = —r 2 dr dcos(8) dip and 6(f— s) = 4j S(r — s) 5(cos9 — cos£) 5(ip — ip). We 
immediately deduce that 

6(r-s) = ±~S(r-s) £ £ [^(CVO]* (*>¥>) (B-13) 

i=0 m=-Z 

On more general grounds it can be shown that the solutions of any Sturm-Liouville spectral 
problem always form an orthonormal and complete set in the background vector space of the 
problem. Thus, the semi-infinite set of the Legendre polynomials {Pq}i=q,i,2,-- are the solutions 
for the eigenvalues A = 1(1 + 1) of the Sturm-Liouville spectral problem 

{d x [(l~x 2 )d x }+A}f(x)=0 (B.14) 

on the real interval [—1,1] with the boundary condition that / must be finite everywhere on 
[—1,1]. As such this set then is orthonormal and complete in the vector space ji of square 
integrable functions on [—1, 1] and allows for series expansion in this space. We more concretely 
have 

nl / i\ , , , A/, i 



-i 



I]'Arll]U.r) [! + -) dx = 5 lk | ^1 + ^ );>',(, -)./>',(// > = ,)(,■-,/) (B.ir,) 



When dealing with the functions /(ri,r 2 ), which depend on two vectors ri(r\, 9±, ip±) and 
^2( r 2j 9 2 , if 2 ), the multipole expansions must be performed over the collection of bi-polar spher- 
ical harmonics. These are obtained by irreducible tensor product of the spherical harmonics 
with different arguments 

(Y h (9 l ^ l )®(Y l -(9 2 , V2 )) LM = £ (hm 1 l 2 m 2 \LM)Y^ 1 (9 1 ,^ 1 )Y^(e 2 , i p 2 ) (B.16) 

mi ,m2 

and form a complete orthonormal set in £| 2xS2 . lJ,\m\l 2 m 2 \LM) is a Clebsch-Gordan coef- 
ficient (cf. Appendix IC.3[) . Generalization to the functions of more than two vectors and to 
multi-polar spherical harmonics is straightforward [17]. If f(r\,r 2 ) is invariant under rotation 
of coordinate systems then it depends only on r\, r 2 and r*i ■ r 2 = r\r 2 COSW12, where coso;i2 = 
cos 9\ cos 9 2 + sin 9\ sin 9 2 cos(y>i — ip 2 ), and its multipole expansion contains only the zero rank 
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bi-polar spherical harmonic (Y^ </Ji)®Y i2 (0 2 , V>2))oo = j 2 tl+i ( Y 'i (^2, ^2)) ft^, 

where the scalar product 

2/ + 1 

Y ; (fli,^i)-Y ; (^ 2 )= J] [y^(fli,^i)] y i {,(fl2,^) = -^-P^(cosw 1 2) (B.17) 

m=—l 

The last equality expresses the addition theorem for the spherical harmonics, proven as follows: 

1 

P(5(cosw 12 )= ^ a* m y^(0 2 ,y 2 ) 

i m— — l 

+ j- J J sin 6 2 d0 2 J dif 2 [Y^(e 2 ,ip 2 )]*Y l (uJi 2 ,ji 2 ) 

but Yl{6 2 ,ip 2 )= &J,y>i2,7i2) r» ¥>i) , = , ^(^V and 

&;o = y sin lu 12 duj 12 J d-y 12 [Yj n (u 12 ,'y 12 )]*Y^(e 2 ,ip 2 ) so a; TO = ( ^ + 1 J [^m^i^i)]* 



the angular variables (wi 2 ,7i 2 ) and {9 2 ,ip 2 ) in the evaluation of bio being exchangeable. 

Let now us focus our attention on the series expansion of exp(ia;y) over the Legendre polyno- 
mials P&(y) : 

OO 

exp(ixy) = Y / ci(x) p o(y) ( B -!8) 
ci(x) = eM™y)Po(y) (l + 5) dy = ^±-^ cxp ( lxy )[d l y (y 2 - 1)'] dy 
Integration 1 times by parts we get 

where 

J v {x) = 1 (^Y [ exp(iaj,)(l - y 2 y-? dy (B.19) 

7T2 (y — -y. \Z/ J_! 

J v (x) defines the Bessel function of index v and is a particular solution of the Bessel equation 



,2 



f(x) = (B.20) 



The general solution Z v (x) can be apprehended on representing the function g{x) = x v fix) 
by the generalized Laplace integral 



g(x) = / cxp(xz) v(z) dz 
Jv 
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and searching for analytic functions v(z) and paths V in the complex plane that verify the 
equation [xd 2 + (2v + l)d x + x] g{x) = 0, that is such that 

= (2v + 1) / exp(xz) zv(z) dz + exp(xz) x(l + z 2 )v(z) dz 
J-p J-p 

= (2v + 1) / exp(xz) zv(z) dz — j exp(xz) d z [(1 + z 2 )v(z)] dz + | exp(xz) (1 + z 2 )v(z) \<p 



or equivalcntly 



zu( z ) 

exp(a;z) [{2v + 1) 2 — ^ — d z u(z)} dz — | exp(xz) u(z) \p> = 



where u(z) — (1 + z 2 )v(z). It suffices that the integral and the integrated term in this last 
expression are separately null. This implies that [2v + l)pr+i ~ d z u{z) — 0, that is u(z) = 
C{l + z 2 Y +1 2. So 

Z v (x) = Cx v / exp(a;z) (l + z 2 ) u ~i dz (B.21) 



with a path V such that the variation of exp(ccz) (1 + z 2 ) ly+ 2 over it is null. The solution 
Z v {x) = Jv(x) is deduced by observing that (1 + z 2 ) v+ ? = for z = ±i provided that 
5i(z/) > — |, which suggests to take for V the straight line from —i to i and to perform the 
z = iy variable change. A recipe to retain is that the method of Laplace integral in the complex 
plane works for all the homogeneous linear differential equations where the coefficients are linear 
functions of the variable. 

x~ v J v {x) is an analytic and even function of x, which thus admits a series expansion every- 
where convergent over even integer powers of x. This allows forming a series for J v {x), which 
necessarily is convergent, and by analytic continuation defining J v (x) for any v. Concretely, 



J ^xp(ixy)(l-y 2 )»~? ^ = 2 ^ 



cos(xy)(i -rr 2 d V 



k=0 

(cf. eq. (|B.3|) and eq. (|B.5|) 1. therefore 



2 2k k\{v + k) 



Let r be a vector in the real space R 3 with spherical coordinates (r,9 r ,tp r ), q a vector in the 
dual space R 3 with spherical coordinates (q,0q,(p q ) and w the angle between them. If qr is 
substituted for x and cos a; for y in the expansion of exp(ixy) over {Pq}i=o,i,2,--- then, using 
the addition theorem for the spherical harmonics (cf. eq. (IB.17j) ). we find 



exp(i^f)=47r5^* , i J (gr) ]T [Y^fa, <p q )} * Y l m (6 r , <p r ) 



(B.23) 
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where ji(x) = ygj Ji+i ( x ) defines the spherical Bessel function of index I. As from the series 
expansion of x~ v J u (x) over even integer powers of x (cf. eq. (|B.22j) ) and by making use of the 
identity (cf. eq. (|B.5[> ) 

V^(2z + l)!-2 2z+1 z! (z + iy 

it is inferred that 

sin a; 2 l V 

jo( x ) = —^r and ilW»->o ~ + \y x ( B - 24 ) 

Also useful is the recursion relation 

j^ix) + Jl+1 (x) = + j,^), Z > 1 (B.25) 



which also is deduced from the series expansion of x v J v (x) over even integer powers of x 
(cf. eq. (|B.22|> ) by merely observing that 



k=0 



k=0 



°° 1 I / T \ 

fc=0 v ; 



= ±(-r^ m „l 1+h)! (f ) ,+2 * - ---ww 



" ¥ (hM(u 

h=0 

so that 



d x {J v {x)) + -J v (x) = J„-i(x) and d x (J v (x)) J v (x) = —J v+1 (x) 

X X 



C ANGULAR MOMENTA 
C.l Definition and Properties 

Angular momenta by definition arc infinitesimal generators of the state transformations asso- 
ciated with the spatial rotations. We shall symbolize these 

|^)^|V') = U(a;,fl)|V') (CI) 

for a quantum system subject to the rotation R(u>, u) about the unit vector u through an angle 
ui. TJ(u>, u) materializes a mapping of the state space of the quantum system over itself which is 
bijective and preserves the modulus of the scalar product so, by virtue of a Wigner theorem, is 
either a unitary or an anti-unitary operator up to phase factors. We can exclude the eventuality 
of the anti-unitarity for this changes the sign of the commutators, which would be inconsistent 
with the infinitesimal transformations, and we can fix the phase factors on imposing that the 
set of the operators XJ(u>,u) equipped with the composition law for the state transformations 
form a group Q homomorphic to the group G of the spatial rotations. We call to mind that the 
kernel K, of the homomorphism is an invariant subgroup of Q and that the quotient group Q /JC 
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is isomorphic to G [1] . Anyway, no phase factor should be expected without contradiction for 
the infinitesimal rotations R(5uj 1 u), so that lim^^o U(w, u) = 1 Vu which in turn implies that 

3 Ju : TJ(6w, u) = 1 - i Slo J fi + O(Slu) 2 (C.2) 
Ja is self-adjoint since U(5w, zt) is unitary. Vv R(Sw, u) v = v + 5u)(u A v) + 0(Suj 2 ) so that 
R(5p,n)R(5u,u) v = v + 5lj(u A v) + 0(5lu 2 ) + 8p(h A v) + O(SpSuj) + O(Sp) 2 
and, Bi being basis vectors, R(5u, u) = R(Slj, J2i(^i ' fyii) = Hi[R((ii ■ u)Su), e*)] + 0(5lu) 2 . 
We deduce that 

J a = ^(e 4 • u)3 &i = it ■ J (C.3) 

by interpreting the operators as the components <5j • J of a vector J over the basis (<5j). J 
defines the total angular momentum of the quantum system under concern. A finite rotation 
can be built up from successive infinitesimal rotations about the same axis, to be precise 
R(uj,u) = limjv^oo R(^,u) N . Accordingly, 

/ 2\ N 

\J(cj,u)=lim ( 1-i ^ («■ J)+o(^) J = exp{— i lu (ii ■ J)} (C.4) 

Notice that U(o> + = U(5u), u)U(u), u) = [1 — i 5co (u • J)]U(cj,m) so <9 w U(w,u) = 

— i (it • J)U(w, m), which, together with U(0, u) = 1, is solved as U(co, it) = exp{— i cj (u ■ J)}. 

A finite rotation is often described as the succession R(uj,u) = R(j, Z)R((3,y')R(a, z) of a 
rotation about the z— axis through an angle a — (y,y'), a rotation about the new y' — axis 
through an angle /? = (z, Z) and a rotation about the new Z— axis through an angle 7 = (j)', Y). 
a, (i, 7 are called the Euler angles of the rotation, which then is denoted R(a, (3, 7). U(w, u) as a 
function of these angles writes U(a, (3, 7) = exp{— i 7 (Z- J)} exp{— i /3 (i/'-J)} exp{— i a (£• J)}. 
The same rotation is obtained by first performing a rotation about the z— axis through an angle 
7, then a rotation about the initial y— axis through an angle /? and finally a rotation about the 
initial z— axis through an angle a, that is we also have R(oj, u) — R(a, z)R{(5,y)R{^f, z) and 

U(a, /?, 7) = cxp{— i a (z ■ J)} exp{— i [3 (y ■ J)} cxp{— i 7 (z ■ J)} (C5) 

Let O be an operator acting on the states \ip) of the quantum system. If this is submitted 
to the rotation R(lu,u) then V <V|0|^> = [U(w, u)] O U(w,m)+|V') = (V'lO'lV)- Wc 
deduce that the operator transformations associated with R(cj,u) are given by 

0^0' = U(w,m) O [U(w,u)]+ (C.6) 

S is a scalar operator if and only if Vw Vu U(w, u) S [U(w, it)] + = S and Vg ( is the component 
ii ■ V of a vector operator V if and only if Vu U(a>, u) Vj, [U(ct>, = (R(w, u)e») • V = 
• V. When the rotation is infinitesimal the operator transformations take the form 

O -> O' = O i (Jw [ u • J, O ] + 0((5w) 2 (C.7) 

We deduce that V is a vector operator if and only if Vw Vw [ u • J, w • V ] = i (u A 6) • V, 
since u-V — > w • V — i <5a; [ u • J, ■ V ] + 0(<5w) 2 but also w • V — > (R(uj 7 u) v) ■ V = 
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(v + Slu(u A v) + 0(5uj 2 )) ■ V. An alternative definition of angular momenta is extracted from 
this, namely as vector operators J satisfying the commutation relations 

[it -3, v- J ] = i (it A if) -3 (C.8) 

where (it, v) is any pair of ordinary vectors or of vector operators that commute with each other 
and with J. Notice that this is a more general concept, which coincide with that of the total 
angular momentum of the quantum system under concern solely if for all the associated scalar 
S and vector V operators we further have [it- J, S ] = and Vu Vu [ u ■ 3, v- V ] = i (it A V) ■ V. 

Using the covariant e m or contravariant e rn — e m Cartesian basis vectors in place of the vectors 
u and v the commutation relations defining J are merely expressed in the form 

[ J m , J„ ] = i e mnl 3i (C.9) 

where 3 m — e m ■ 3 (m — x, y, z) are the Cartesian components of J and where e mnl = 1 if 
(m, n, I) is co-cyclic to (x,y,z), = —1 if (m,n,l) is anti-cyclic to (x, y, z) and = otherwise. 
Using the covariant spherical basis vectors e_i = ^(e^ — i e y ), io = e z , e+i = ~^(^x + i e y ) 
we get 

[ J MJ 3 U } = -y/2 (1/*1i/|1A) Ja (Ai,i/,A = ±1,0) (CIO) 



where 



3-i= e-i-3 = -j=(3 x -i 3 y ), J =e -J = J 2 , J+i = e+i ■ 3 = --j=(3 x + i 3 y ) 



are the covariant spherical components of J and a Clebsch-Gordan coefficient (cf. Appendix lC.3j) . 
Using the contravariant spherical basis vectors e M = e* = (— ) A 'e_ AI we rather get 

[ J M , 3 V ] = V2 (l/xli/|lA) J A (/*,!/, A = ±1,0) (C.ll) 

where J M = [J M ] + = (— ) /i J_ A1 are the contravariant spherical components of J. We shall take 
this opportunity to call to mind that e M • e v = <5 M- „ (= 1 if /1 — v and = if /1 ^ v) and that 

if A = E M = ^e" and B = £ M = £ p B M e" 

then I- B = = E M (-)"^ B -P 

whereas {A A B}\ = -i \/2^2 ^ v {l^lv[l\) A^B U and {A A B} x = i ^(l/ili/llA^B". 

Whatever the basis vectors ej it is effortlessly inferred that [J 2 , e-i ■ 3 ] = 0, in particular 

[J 2 , e M -J]=0 (/i = -l,0,l) I [ J 2 , e TO • J ] = (m = x,y,z) (C.12) 

J 2 commuting with Jo share with it common eigenstates \tJM), which, because the operators 
are self-adjoint, are orthonormal : (tJM[t' J' M') = 8 t ,t'8j,J'&m,m> ■ Using the commutation 
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relations [ J M , 3 U ] = — \/2 (l/xli/|lA) J a and [ J 2 , ] = (/i, v, A — ±1, 0) together with the 
fact J 2 = X)u( — ft i s an eas Y matter to show that (cf. the case J = Lin Appendix IB]) . 



3 2 \tJM) = J(J+l)\rJM) 

J \ T JM) = M\tJM) 
3 +1 \tJM) = -X m \tJM + 1) 
J_i|rJAf) = X M -i\tJM- 1) 



(C.13) 



where X M = [\{J - M)(J + M+ 1)] 2 with 2J e N and 2M G Z : —J < M < J, that is 
J = 0, |, 1, |, 2, ■ • • and, J being fixed, M = — J, — J+ 1, • ■ ■ , J— 1, J. r distinguishes between 

the orthogonal eigenstates of J 2 with the same eigenvalue J(J + 1). The set {\tJM)}-j<.m<j 
engenders a state subspace iS T ,j of dimension 2 J + 1, which then admits the closure relation 



J 



J2 \tJM){tJM\=1 tJ (C.14) 



M=-J 



It finally may be shown using the same method as for the spherical harmonics (cf. Appendix [B| 
that the set {\tJM)} t jm is complete in the state space S of the quantum system under concern. 



C.2 Wigner D-Matrix 

The matrix representatives of the operators {J^}^=-i,o,+i over the basis {\t JM)} t jm, trace- 
less and hermitian, materialize the irreducible matrix representations of the su(2) Lie alge- 
bra. These ascend faithfully to the associated SU(2) Lie group, because this is simply con- 
nected, namely the irreducible components S T .j of the representation space S are also those 
of SU(2) and the matrix representatives of the operators XJ(uj,u) or U(a, /?, 7) over the basis 
{\tJM)} t jm provides with the irreducible matrix representations of SU{2). Using the Euler 
angles 

(t'J'M'\ U(a,/3, 7 ) \tJM) = (C.15) 
= (t'J'M'\ exp{-i a (z ■ J)} exp{-i (y ■ J)} exp{-i 7 (z ■ J)} \tJM) = 
= $t',tSj>,j exp{-i aM'} d J M , M ((3) exp{-i 7M} = 6 T ^ T 5j,j T> J M , M {a,(3,^) 

T)f IIM (a, (3, 7) is sometimes dubbed a Wigner D-function and T> J (a, /?, 7) a Wigner D-matrix. 



J = |<r for J = h, where <x is the Pauli operator. The matrix representatives of its Cartesian 
components, (x ■ <r) , (y ■ <x) , (z ■ er), over the basis {\t^ + 5), |r| — i)} are the Pauli matrices 




(C.16) 



(T 2 == <jy = fx 2 = 12x2 (unit 2x2 matrix) and <J m a n = —cr n a m = iai for any cyclic permutation 
(to, n, I) of (x, y, z), which allows showing that 

(t? • a)(v ■ d) = (u ■ v) 1 + i <r ■ (u A v) (C.17) 

for any pair (u, v) of ordinary vectors or of vector operators commuting with <r but not neces- 
sarily with each other. In particular, Vu (u ■ er) 2 = (u ■ u) 1 = 1 (u : unit vector). So 

exp{-i ^ {u ■ <?)} = cos(-) 1 - i sin(-) (u ■ <r) (C.18) 
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and 



Z>*(a,/3, 7 ) = 



cos(|) exp{-z ^} - Biii(|) ex P{"* ^> 



sin(f)exp{z ^p} 



cos(f)exp{« ^} 



(C.19) 



V J (a, /3, 7) or rather d J M , M (f3) for any J is elegantly evaluated from X>2 (a, /3, 7) using the spinor 
formalism, which, for each fixed r, allows building up the angular momentum states \tJM) 
from tensor products of the spinor basis states x± = exp{i (^±}|ri ±i). As to avoid writing 
heaviness the tensor products will be denoted multiplicatively, for instance Xl" ' XN — IL=i Xl 
for O7+X+ + ViX-) ® ••• ® (?7nX+ + %-X-) = ®jIi(^ + X+ + CX-)> ^ G C, and x° for 
X ■ ■ • X- One of the essential tools of the spinor calculus is provided by the mappings 

JV 

d x± ■ Xl • ■ ■ XN d x± (xi • • ■ Xjv) = X! ' ' ' Xl-i (X± I Xi) Xl+i -"Xn (C.20) 

1=1 

which might be interpreted as derivation operators. On examining the action of the components 
of the vector operator J on the states x± = exp{i y±}|Ti ± |) (cf. eq. (|C. 13|) ) the following 
identities are inferred 

J-i = ^ ex P{* ~ V-Wx-d x+ , Jo = 7^(x+d x+ -x-d x _), (C.21) 



1 

1 

2 ( 



J + l = 7r cx P{-« [<£+ - ^p-WX+dx- 



J 2 =K(K + 1) with K = -{ X +d x+ + X -d x _ 



Using these it easily is shown that Va e N V6 e N Jo(x+X-) = 5( a _ &)(X+X-) and 



J 2 (x+X-) = { 9 4^-){ 9 4^- + l)(x+X-)i that is any monomial x+X- is an eigenstate of J and of 



J 2 simultaneously. It also is found out that J-i(x+X-) = -^ acx p{* [<P+ — ( P-]}(x+ X X- ) 
and J+i(x+X-) = -7f &ex P{-* [<P+ ~ V-lKx+^X-" 1 ), that is applying J_i and J +1 on 
any monomial x+X- allows generating the irreducible component <S a+t of the representation 
space S. All this merely suggests that 

\tJM) = c(j,m) x j + +m x j -- m 

If C(J, J) is fixed to -j=L= then, since [J^^tJJ) = (^) n ^0^\tJJ - N), we find 

(e^ifa -*-]})'-» 

y/(J + M)\(J-M)\ 

by applying J — M times J 1 on \tJJ). We are free to also fix the phase (<p + ~ <p_) with 
some arbitrariness. A conventional choice is (ip+ — <ys_) = (2rn + with m € Z, so that 
exp{i [(p + — <p-}} = —1 and 

J+M J-M 

WJM) = l-) J - M — ^+ x - (C.22) 
V(J + M)!(J-M)! 



Prl-36 



JOURNAL DE PHYSIQUE IV 



U(0, /3, 0)|r JM) now can be computed either in the form 

- j-m (X+ cos I - X- sin |)' 7+M (X+ sin | + X- cos f ) 



(-)- 



§_\J-M 

2 - 



V(J + Af)!(J-Af)! 



(C.23) 



£ |rJM / )(rJM'|U(0, J 8,0)|rJAf)= ]T (-) J ' 



M'=-J 

Using the identities 



M'=-J 



J+W J-M' 

-M A+ X- 

^(J + Af')!(J-M')! 



4mW 
(C.24) 



^(i±^ = (±r 7 -i— (i±xr 



k=0 



k\{n-k)\ 



x k y n ~ k 



(n-l)\ 

n 

and d2(fg) = J2 



k=0 



k\(n-k)\ 



d k f d n x ~ k g 



we deduce 



^m'm(^) — ( )' ^7 



(J-M')! 



M' + M Af ' -jVT 

(l + cos/3) 2 (l-cos/3) 2 x 



{J + M')\{J + M)\(J - M)\ 

xdLtf [(l + cos/3)^ M (l-cos/3) 



= (-) ^ 



(J + M)! 



M' + M , , M'-M 

(1 + cos^) 2 (l-cos/9) 2 x 



_(J-M)!(J + M')!(J-M')! 

x [(! + cos/3) J+M '(l - cos/3) J - M '] (C.25) 

then, from this, the expression of T> J (a, f3, 7) for any J. We find the symmetry properties 
d M >M{-P) = d J MM >(0), d J _ M ,_ M {p) = (-) M '- M d J M , M (0), 
d J MM ' W = (-) M '- M <'m(/3), d J _ M _ M , ((3) = d J M , M (/3) 

<, M (/3±2mr) = (-) 2nJ d J M , M (p) (n G N), 
<, M (/?± (2n + 1)tt) - (-) ±(2 " +1)j - m <^m(/?) (n G N) 



Setting J = / G N, M' = m G N (-/ < m < I) and M = Q or J = Z G N, A/' = and 
M = m G N (— Z < to < Z) in the eq. (|C.25[) and comparing with the functional forms of the 
Legendre polynomials P l m and of the spherical harmonics (cf. eq. (|B.7p ) it is found out that 

dlo(P) = R m y!^| ^i(cos/3) I d l Om (0) = ^(cos/3) (C.26) 

so that 

^o(«,A 7 ) = ^Ill(A-«) I V l 0m (a, 0, 7) = ^ ^ Y^-ft -7) (C27) 
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If R((X3, /?3,73) = R(ct2, /?2 5 72) , /3i , 7i ) describe two consecutive rotations then the Wigner 
function V J M , M (a s , fo, 73) = J2m»=-j v m>m" ( a 2, #2, 72) U^/ M ,,(ai,/3i,7i) whereas the Eu- 
ler angles as, ^3,73 are given as functions of the Euler angles ai,/3i,7i and a 2 , ^2,72 by 
cot(a3 — a 2 ) = f{/hifli)> cos/?3 = cos/?iCOS/3 2 — sin/3i sin/3 2 cos(ai +72) and cot(73 —71) — 
f(0i,/3 2 ), where f(x,y) = cos x cot(ai + 72) + cot y sin( s "^ 72) , whence by setting J = !eN and 
M' = M = we recover the addition theorem for the spherical harmonics (cf. eq. (|B.17jl ). 

The rotation R(tv,u) can be described as the succession R(a^ = tpai (3% = —0u,"/3 = —<Pa) 
R(ot2,/32 = 0,72 = 0) R(oti = (fu, Pi = #a,7i = — (fa) of a rotation (1) which align the 
z— axis along u(9a, fa), then of a rotation (2) about u(9a, tpa) through an angle lo and finally 
of a rotation (3) which send back the z— axis to its initial direction. Using these rotation angle 
uo— rotation axis u(9a, (pa) variables (t'J'M'\ \J(u>, u) \t,JM) = <5 T ' :T <5 U J M , M (yi, 9a, <fu) with 

J 

U J WM {tt,9a, Va) = ^2 ^M' M "( l Pu,0u,-ipu) exp{-iM"u>} V J M , IM {ipu,-0u,-iPu) (C.28) 

M"=-J 

Another way obtaining U'I iim {lo, 9a, (pa) is by the direct (a, (3, 7) — * (cu, 9a, (fa) variable substi- 
tution in T> J M , M {a, f3, 7), with the aid of the relations 



.(I) = cos (§) cos (£±z) , tan(0«)=tau(f) [am(2±X)] 1 



— 2 ^ 2 



M'+M 



^, M (w, ^) = i M '~ M exp{-i(M' - iW> fi } I i itan2C ° = I 



'1 + tan 2 f cos 2 1 



where £ is determined by sin | = sin ^ sin 9a- 



(C.29) 



The trace Tr \U(w, u)] of the operators U(w, u) = exp{— i w (u- J)} over the basis {\tJM)} t jm 
provides with the characteristic functions, or simply, the characters of the SU(2) irreducible 
representations. Since Tr(Pc(Ai ■ ■ ■ An)) — Tr{A\ ■ ■ ■ An) for any cyclic permutation Pq of 
the operators A\ • • ■ An, the characters are invariant under any unitary transformation. As a 
result these are independent of the colatitude 9^ and longitude (pa of the rotation axis, to be 
precise Tr [U(w, u)] = x J (w) with 



x »= £ e Xp {-iMu} = Sin = ^rrflUt cob [(2J + l)g 



(C.30) 



U being unitary Tr [U = Tr [U(a;,w)]. x J ( w ) is a rea l ([x' 7 ( aj )]* = X' 7 ( w )) an d even 

(x J (— = x' 1 ^)) function of u, which is 2tt— periodic if 2J is even and 2n— antiperiodic if 2J 
is odd (x J (w + 2tt) = (-) 2J x J M). 

The characters x J (w) of the irreducible representations of SU(2) can be generalized in the 
form 

X» = V£T+l ^^|^ (ain|)' V(«) (C.31) 

where A € N, < A < 2 J. Xa( w ) ^ s dubbed the associated character of order A of the 
irreducible representation of rank J. xi{ w ) = X J ( W ) an d [xa( w )]* = Xa( cj ) = ( — ) A Xa(~ w )- 
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The collection {Xa( w )}./>^ ^e[o 2-n] exhibits the orthonormality and completeness properties 

f\in^ X t{u) X {>{u)<L>=v8 JuJ2 | £ XaV^M = (C32) 
Jo 2 ^ sm -f 

Using the generalized characters the function U^ iim (lu, 9u,tpu) may be expanded in a series of 
spherical harmonics depending on the angles (9^, ip^) of the rotation axis : 



U&, M (u, 0«, <p*) = EH) A ^7 XaH (JM'MJM) ^/^y Y*(du, Vu) (C.33) 

A, /a 

It then is seen that U^, m (oj, 9-a,ipu) depends on M' and M only through the Clebsch-Gordan 
coefficients (JM'Xfi\JM) (cf. Appendix [03]) [17]. 



C.3 Clebsch-Gordan and 3jm Symbols 



Let Si and iS 2 be the state spaces of two kinematically independent quantum systems or parts 
of a quantum system. The state space of their union is the tensor product S — Si(£)S2- 
Under the spatial rotation R(u>, u), the states \tpk) S Sk (k = 1, 2) get transformed into \ip' k ) — 
Ufe(o;, u)\ipk) = exp{— i a; (u • Jfc)}|^fc), where Jfe (fe = 1, 2) defines the total angular momenta 
of each of the two quantum (sub-)systems, therefore the states <£> \ip2) € <S get transformed 
into \ip[)®\ip' 2 ) = Ui(cj, ®U 2 (w,fi)|^) = (Ui(w,fl) O li)(l 2 <8> Ui(w,fi))|Vi)<8)|^) = 
(Ui(o;,u) <g> U 2 (<J, u))|^/>i) <g> 1^2), but the whole quantum system is featured by its own total 
angular momentum J so that \ip[) ® \tp 2 ) = U(w, ® IV^) = exp{— i w (u • J)}|^i) eg) |i/> 2 ). 

As a result, U(u;, ft) = Ui(w,ft) £g> U 2 (u;, ft), namely 



exp{— i a; (ft • J)} = exp{— i u> (ft ■ Ji)} (B> exp{— i a> (ft • J 2 )} 
and in the case of an infinitesimal rotation R(Slo, ft) 



(C.34) 



1 — i Slu ft ■ J + 0(<5a;) 



so that 



i Suj ft ■ Ji + 0((5w)^ 



Ji 



1 2 — i 5uj u ■ J 2 + 0(8u>y 



(C.35) 



Usually denoted J = Ji + J 2 par abus d'ecriture, this defines the addition of two angular 
momenta. Owing to the kinematic independence of the two quantum (sub-)systems, [ Ji , J 2 ] = 
0. As an outcome, J shows the commutation relations of an angular momentum. Additionally, 
[ (Jfe) 2 , J 2 ] = and [ (Jfe) 2 , Jo ] = (k = 1, 2), which implies that the eigenstates common to 
J 2 and J (cf. eq. (|C.13[0 . engendering the irreducible components S Tt j of the representation 
space S, are eigenstates of (Jfe) 2 with the eigenvalues Jk(Jk + 1) (k = 1,2), hence denoted 
\tJ\ J2JM) (t = T\T2f])- These form a complete set in S. The eigentates |rfe JfeMfe) common to 
(Jfe) 2 and (Jfe)o, engendering the irreducible components S Tk j k of the representation space <Sfe, 
also form for each k a complete set in Sk (k — 1, 2), so the states \t\ J\Mi) <g> |r 2 J 2 M 2 ), which 
we shall denote |nr 2 J\M\ J 2 M 2 ), form another complete set in S = Si (££>iS> 2 . J 2 is computed 
equal to (Ji) 2 + (J 2 ) 2 + 2 £ u (-)"(Ji) M (J 2 ) M , so that [ J 2 , (J fc )o ] + {k = 1,2), whereas 



Jo = (Ji 



(J 2 )o trivially, so that |tit 2 J\M\ J 2 M 2 ) is an eigenstate of Jq with the eigenvalue 



M = M 1 +M 2 



(C.36) 
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Let N TlT2 (J) be the number of sets {\tJ\ J 2 JM)}_ j<m< j associated to each eigenvalue J of 
J 2 and let n TlT2 (M) be the degeneracy attached to each eigenvalue M of J , both for fixed 
quantum numbers tit 2 . n TlT2 (M) — J2.j>\m\ N T it 2 {J) then N TlT2 (J) = n TlT2 (J) — n TlT2 (J + l), 
but n TlT2 (M) is merely the number of pairs (Mi, M2) such that M — M1+M2. As to enumerate 
them it is helpful to take them as points with the abscissa X — M\ and ordinate Y — M 2 
in the plane. n TlT2 (M) then is the number of points in the diagonal X + Y = M. We find 
n TlT2 (M) = if |M| > ,h + J 2 , n TlT2 (M) = J x + J 2 + 1 - \M\ if Ji + J 2 > |M| > \.h - J 2 \ 
and n TlT2 (M) = 2min(Ji, J 2 ) + 1 if |Ji - J 2 | > |M| > 0, which implies that N(J) = 1 if 
J = Ji + J 2 , Ji + J 2 — 1, • • • , I Ji — J 2 | and iV(J) = otherwise: the quantum number r\ is 
limited to a single value then useless (r = tit 2 ) while J must satisfy the triangular condition 

\Ji - J 2 | < J < J\ + Ji (C37) 

Let U Jk (w, w) be the matrix representatives of the transformation operators Ufc(w, u) over the 
basis {\TkJkMk)}-.j k <M k <j k of S Tk .j k , defining the irreducible matrix representations of SU{2), 
for each quantum (sub-)systems (k — 1,2). The tensor product U Jl (u),u) ®U ,,2 (u 1 u) defines 
a matrix representation of SU(2) for the whole quantum system. Its character is computed to 
be 

Tr [U Jl (uj, u) ® U' h (w, fl)] - Tr [U' h (w, u)] Tr [« 72 (w, u)] = 
Ji J 2 J1+J2 J 

= exp{-i Mi lu} x cxp{— i M 2 uj} = J3 ^ cxp{— i M : 

Mi=-Ji M 2 =-J 2 J=\Ji-J 2 \ M=-J 

that is 

X Jl (^)X J2 M = J i J 2 J > (C38) 
j 

where { Ji J 2 J} = 1 if | Ji — J 2 | < J < Ji + J 2 and {Ji J 2 •/} = otherwise. As a corollary, 
the sum of characters of ranks Ji, J\ + 1, Ji + 2, • • • , J 2 may be factorized according to 

53 x J M = x^Mx^M (c.39) 
j=ji,ji+i,... 

U Jl (u),u) ®U j2 (u),u) — ® J j=\jI_j 2 \ U J {ui, u), where U j (lj,u) designate the matrix represen- 
tatives over the basis {|t Ji J 2 JM)}_j<m< j of the transformation operators U(w,u) defining 
the irreducible matrix representations of ST/ (2) associated with the whole quantum system, 
and S TliJl ® 5 T2;J2 = ©£}£_ Ja | 5 T)J . Incidentally, (2Ji + 1)(2J 2 + 1) = ££}£_ Ja |(2.7 + 1) 
which is easy to directly confirm. 

The sets {\tJ\J 2 JM)} and {|r J\M\ J 2 M 2 )}, orthonormal and complete in S, are related to 
each other by a unitary transformation: 

Ji J 2 

\rJxJiJM) = J2 (\tJ 1 M 1 J 2 M 2 )(t,J 1 M 1 J 2 M 2 \) \t.J x J 2 JM) 

\tJiM 1 J 2 M 2 )= (\tJiJ 2 JM}(tJiJ 2 JM\) \tJiM!J 2 M 2 } 

J=\ J1-J2I M=-J 

where the {t J\M\J 2 M 2 \t J\J 2 J M) coefficients arc dubbed Clebsch-Gordan. Complex con- 
jugate to (j J\J 2 J M\t J\M\J 2 M 2 ) they manifestly are r— independent, then rather denoted 
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( J\M\ J 2 M 2 | JM). These can be evaluated in a variety of ways exploiting unitary prop- 
erties and recursion formulas or by operator techniques based on commutations or else by 
means of group theoretical methods. Use can also be made of the spinor formalism, taking 
advantage of its formal analogy with functional analysis to first determine ( J\M\ J 2 M 2 | J J) 
then deduce ( J1M1 J 2 M 2 | JM) from the application of J^7 M on both side of the equation 
\t J1J2JJ) = Y.m u m 2 ( J i M i J ^WJ) \tJiMiJ 2 M 2 ): let ( X k)± = ±|t| ± |) be the spinor 
basis states and dt Xh )± the associated spinor derivative operators of each of the two quantum 
(sub-)systems (k — 1,2) for fixed quantum numbers t\t%. It is inferred from 

J+i|r J1J2JJ) = [(Xi)+0( X1 )_ + (X2)+d ix2) _] \tJiJ 2 JJ) = (C.41) 

(cf. eq. (|C.21[) with (<pk)+ — ((fk)- = and from our knowledge of homogeneous partial 
differential equations that 

\rJ.J2JJ) = f ([(xi)-(xa)+) - (Xa)-(xi)+] , (xi)+, (xa)+) (C42) 

Additionally, \tJiJ 2 JJ) is an eigenstate of J| = K/ £ (Kfc + 1) with eigenvalues Jfe(J/c + 1) thus 
also of Kfc with eigenvalues We then get an Euler partial differential equation 

(Xk)+d (Xk)+ + (Xk)-d ixk) _f = 2J k f (C.43) 

(cf. eq. l|C.2ip ). This implies that / is an homogeneous function of order 2Jj~ in (xfc)±> that is 
of the form 

\rJ,J 2 JJ) = £ C L [(xi)-(Xa)+) - (x 2 )-(xi)+] L (x^ 1 "^) 2 / 3 ^ (C.44) 

L 

It follows that Y.l -J1-J2 + J)[(xi)-(X2)+) - (X2)-(xi)+] L (xi)+ Jl " L (X2) 2 + J2 - L = 0, 

since 



MTJ1J2JJ) = l\J2 

I fc=l,2 



(Xfc)+ - ( Xk )-d (Xk) _] I |rJi J 2 JJ) = J\rJxJ 2 JJ) (C.45) 

(cf. eq. (|C.2ip ). Gl then is null unless L — Ji + J 2 — J. Accordingly, 

\rJ,J 2 JJ) = C [(xi)-(X2) + ) - (X2)-(xi) + ] Jl+J2 - 7 (xi)i + ' 7l " J2 (X2)i +J2 - 71 = (C.46) 

= c Jl+ f; J ( Jl + J J- J )\ (xi) x (xi) 2 + /i - x (-) Ji+j2 - / - x (x2) 72 - /+Ji - x (x2)i 2+J - 7i - x 



x=o 



Substituting Ji — Mi for X and Mi + M 2 for J, the summation extends over all the pairs 
(Mi, M 2 ) with fixed J = M x + M 2 (notice that we do have n TlT2 (M = J) = J\ + J 2 + 1 - J). 
Then 

\tJiJ 2 JJ) = C {J\ + J 2 - J)\ x (C.47) 



\- ^ I \Ji-Mx x (J!+Mi)!(J 2 + M 2 )! ,_ T 

L 2^ (-) Wm 2 ,W (Ji _ Mi)!(j2 _ m V-hM^h) 



x 

Mi=-Ji M 2 = -,/ 2 
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since (-) J *- M HXk) J - h - Mk (Xk)l k +M " = y/(Jk - M k )l(J k + M k )\\r k J k M k ) (k = 1,2) (cf. eq. 
Using the binomial identity 



E 



n\ I m 

p) \ r - p. 



n + m 
r 



where 



kj p > kl(p-k)\ 



and 



p<0 



1 ; k\(-p-l)\ [ ' 



k — p — 1 
k 



obtained by considering the binomial coefficients of x r y m+n r in the binomial expansion of 
both side of the equality (x + y) n (x + y) m = (x + y) n+m , it is a matter of elementary algebra, 
putting n = J1 — J2 — J—I < 0, m = J 2 -J\-J-l < 0, p = J1-M1 > and r = J1 + J2-J > 0, 
to prove that 



E 



lM 1 + M 2 ,J 



(Ji + Mi)!(J 2 + M 2 )\ (J + Ji - J 2 )!(J + J 2 - Ji)\(J + Ji + J 2 + 1)1 



Mi=-J u - ,Ji 
M 2 =-J 2 ,-' ,^2 



so with 



(Ji - Mi)!(J 2 - M 2 )\ 



(2 J + l)!(Ji + J 2 - J)! 



(C.48) 



\C\ 



(2J+1)! 



(J + Ji - J 2 )!(J + J 2 - Ji)!(J + Ji + J 2 + !)!( Ji + J 2 - J)! 



|tJi J 2 J J) is normalized to unity. Adopting the conventional choice of phase [12], C is fixed to 
|C|. Now to find the general Clebsch-Gordan coefficients it suffices to operate J — M times on 
\tJ 1 J 2 JJ) with either J_i = ^[(xi)-%. l)+ + (X2)-3( X2 )+] (notice that J_i[(xi)-(x 2 )+) - 
(x 2 )-(xi)+] = 0, which makes the computation easier) or directly with J 1 = (Ji)_i + (J 2 )_i, 
calling to mind that 



(J-x^lrJM) = (-L) w 



(J + M)!(J - M + N)l 
(J — M)\(J + M — N)\ 



\tJM - N) 



We find 



(J-O^lrJ! J 2 JJ) = (4=) J " M 



1 



(2J)!(J-M)! 
(0)!(J + M)! 



\T.hJ 2 JM) 



C (Ji + J 2 -J)! 2 J] ( _)./ 1 -m 1 ,5 Mi+M2iJ1 



M 1 = -J 1 M 2 =-J 2 
J-M 

X 



f (Ji + Mx)!(J 2 + M 2 )! 
(Ji-Mi)!(Ja-M 2 )! 



E 



J — M 
X 



[(Ji)-if [(J 2 )-i] J - M " X |r J 1 M 1 J 2 M 2 ) 



= (j=y-™c { j 1+ j 2 -jy. £ £ (-)■— w 2 ,,-! Jl+Ml)!(j2+M2)! 

^ Mi=-Ji M 2 =-.J 2 



E 



(J-M)! 



X=0,- ,J-M 
Mi-X>-Ji 
M 2 -J+M+X>-J 2 



X!(J-M-X)! 



(J x - M x )!(J a - M 2 )! 
(Ji - Mi + X)!(J 2 - M 2 + J - M - X)!" ' 



(Ji + Mi - X)!(J 2 + M 2 - J + M + X)\_ 

x \tJi{Mx - X)J 2 (M 2 -J + M + X)) (C.49) 
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It finally is inferred, after the substitutions Mi = X + \x\ and M 2 = J — M — X + jj, 2 , that 



(JiHiJ 2 H 2 \JM) = (-) Jl -^ 
"(Ji-Mi)!(^2-M2)! 



(2J + l)(Ji + J 2 - J)!(J + M)!(J - M)! 



(Jl+^l)!(J 2 +M2)! 



(J + Ji- J 2 )!(J + J 2 - Ji)!(J + Ji + J 2 + 1)! 

(J 1 + t x 1 +X)\(J 2 + J-Mi-^)! 
X!(J - M - X)!(Jr - mi - X)!(J 2 — J + mi + X)! 

(C.50) 



with the summation over X such that the arguments of the denominator are non negative. 
Interestingly this is limited to the single term X = when mi — J\ and to the single term 
X = J\ — [i\ when J3 = J\ + J 2 and [i 3 = fii + M2, which provides with the simpler formulas : 



(J1J1J2M- Jx\JM) = 

(2J + l)(2Ji)!(- Ji + J 2 + J)!(Jr + J 2 - M)!(J + M)! 



(Ji + J 2 - J)!(Ji - J 2 + J)!(Ji + J 2 + J + !)!(- Ji + J 2 + M)!(J - M)\ 



(C.51) 



(Jl/(il J 2 /^i 2 | Jl + J 2 yUl + M2) 



(2Ji)!(2J 2 )!(Ji + J 2 + /ii + M2)!(Ji + J 2 - Mi - M2)! 



(2 Ji + 2J 2 )!(J 1 - mi)!(^i + Mi)!(^2 - M2)!(^2 + M2)! 



(C.52) 



As evident from its algebraic expression, i/ie Clebsch-Gordan coefficient ( Ji/ii J2M2I JM) is real. 
It also exhibit a high degree of symmetry. With the substitution X = J — M — Y , for instance, 
we get the result 

(J1MK/2M2IJM) = (-) Jl+j2 - J (J2M2^iMi^M) (C.53) 

A number of other symmetry relations can be found out. These are better displayed by defining 
the 3jm symbol 



Jl J2 J3 

mi m 2 77J3 



= ^_^'3+»"3+2jl 



[k *] (h - m i32 - m 2 \j 3 m 3 ) 



(C.54) 



where [jl 1 ,j 2 2 , • • •] is a standard abbreviation for {(2ji + l) ei (2j 2 + l) 62 • • • } 2 . The inverse 
relation is 

ji ji h 
mi m 2 — m 3 



(jimij 2 m 2 \j 3 m 3 ) = ^ 



The 3jm symbol are null unless the triangular condition \ji — j 2 \ < j 3 < ji + j 2 and the 
zero-sum condition mi + m 2 + m 3 = are satisfied. 

The 3jm symbols are invariant by cyclic permutations of the columns and multiplied by the 
phase factor (— pi+^+Ja by exchange of two columns or by sign change of all the projection 
quantum numbers : 



Jl J2 J3 

mi m 2 m 3 



J3 Ji .12 
m 3 mi m 2 



J2 J3 Ji 
m 2 m 3 mi 



Ji J2 J3 
mi m 2 m 3 



(-) 



31+32+33 



J2 Jl J 3 

m 2 mi m 3 



(C.55) 
(C.56) 
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31 J2 33 \ _ ( _s h+h+h I Ji 32 33 



(-)n+n+33 jl J* J* ( C57 ) 
mi m 2 m 3 I \ —mi — m 2 —m 3 1 



The orthonormal properties of the quantum states 

(jimxjim^hm'Jim'z) = (jimi32 m 2\(J2j 3 ,m ? \33m 3 )(j 3 m 3 \)\j 1 m' 1 j 2 m' 2 ) =^ OTl , m ' l S m2 . m > 2 and 

0'lj2j3"^3|jlj2j3 m 3> = (3l3233m3\(Y, mi ,m 2 32™>2) 0'l TO l J2™2 I) I jl h i 3 TO 3> =^3.^ <Wm£ 

impose on the 3jm symbols the orthonormal conditions 

E ( m \ m 2 TO 3 3 ) ( m \ TO 3 3 ) =*«.i,mi S m2 , m > 2 (C.58) 

e ( & i & ) ( & i 5 a )= w **« «■ * « < c - 59 ) 

where {ji j% j 3 } = 1 if 31,32,33 satisfy the triangular condition and {ji 32 33} = otherwise. 
Since (ji)o + (12)0 — (j3)o is diagonal in both the |T.7imij 2 m 2 ) and \Tj1j2j3m3) representations 
the orthonormal properties of the quantum states materialize themselves for fixed values of 
7713 = —mi — m 2 = m,g = — m\ — m' 2 as well, which provides with the additional formulas 

ym( 31 32 33 )( h , , 32 33 )=s mim > (ceo) 

^—^ \ mi -mi -ma m 3 J \ m\ -m\-m 3 m 3 J 1,1 v ' 

3i 32 h \ ( h 32 3 3 \ _ r -2] 

"ia,j' 3 \JL Jl J3S 

(CM) 



^— ^ I mi —mi — m 3 m 3 ) \ mi —mi — m 3 m 3 / L 3 J J3 J ' 3 ^ -* 



The Clebsch-Gordan coefficients or the 3jm symbols are the basic quantities of angular mo- 
mentum coupling and get ubiquitously into the formulation of a number of physical quantities. 
Let, for instance, V J M , M (ot, (3, 7) = (JM'\ U(a, (3, 7) | JM) be the Wigner D- function of a quan- 
tum system associated with the rotation R(a,(3, 7) (cf. eq. jCT5|) and let m i a i^il) = 

k w 

(JM'\ Ufc(a,/3, 7) |JM) (fc = 1,2) be the Wigner D-functions of two kincmatically indepen- 
dent components of this system. Calling to mind that U(a, /3, 7) = Ui(a, /3, 7) ® U 2 (a,/3, 7) 
and making use of the appropriate closure relations, it effortlessly is shown that 

V J M , M (a,0, 7 ) = £ (JM'\JiM[J 2 M 2 ) 2>£, Mi (a,/3 >7 ) 2^, Ma (a, /?, 7) (JiMiJ 2 M 2 \JM) 

M' lt M 2 
M lt M 2 

^m { mM^^) V m 2 mM^)= E (^M(J 2 M^|JM')^M(a,/3,7) (JM| J x Mi J 2 M 2 ) 



J.M'.M 



Setting J = L € N, Ji = L x e N, J 2 = £ 2 € N and M{ = M2 = in the second equation 
interestingly provides with the frequently used multipole expansion of the product of two 
spherical harmonics. We indeed get 

P \(a,/3,7) 2^>,/3,7) = ^(£i0L 2 0|L0) V% M (a, f3, 7) (LM\LiMiL 2 M 2 ) 

L,M 
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that is 



L,M 



(2L 1 + 1)(2£ 2 + 1)(2L + 1) 

47T 



L\ L 2 L \ [ V Lm\]* I Li L 2 L 



(C.62) 



/ I Mi M 2 M 



(cf. eq. (102711 ) where ft = (-/3,-a) and [Y^(-(3, -a)]* = {-) M Y± M {-(3,-a). Using the 
orthonormality property of the spherical harmonics (cf. eq. (|B.10p ) it also is found out that 



(2£i + l)(2L 2 + l)(2L+l) 

4n > (C.63) 

' L x L 2 L\f L x L 2 L 

J I M x M 2 M 



Specialized formulas for the 3jm symbols can be derived from those obtained for the Clebsh- 
Gordan coefficients. Using the eq. (|C.50j) we compute, for instance, 



(lOjmljm) = , 



whence 



1 3 3 

— m m 



(-) 



j+m— 1 



[r 1 ] 



vW+i) 



(C.64) 



(C.65) 



As from the eq. (|C.5ip we deduce 



31 32 
h ~3i + m 3 



33 

-m 3 



) = (-)"*+*-"* [jg- 1 ] (jijih -h + m 3 \hm 3 ) = 
(2ji)!(— ii +32 + h)Kh + h - m 3 )l(j 3 + m 3 )\ 



. (ji + 32 - 33)1(31 - 32 + J3)!(ii + 32 + 33 + + h + m 3 )\(j 3 - m 3 )\ 

(C.66) 

and from the eq. (|C.52[) 



ji h 3i + 32 
m\ m 2 —mi — m 2 



= [( n + ]2 )- 1 } (j 1 m 1 j 2 m 2 \h+j 2 m 1 +m 2 ) = 



(-) 



-ji+j2—mi—m 2 



(2ji)!(2j 2 )!(ji +32 + mi + m 2 )\(ji + j 2 - mi - m 2 )\ 



_(2ji + 2j 2 + l)!(ji - mi)!(ji + mi)!(j2 - m 2 )!(j 2 + m 2 )\ 



(C.67) 



A few particular instances of these formulas are 

(0 -i =H ' +m[rl] (C - 68) 
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1 3 ~ 1 3 

1 — to — 1 m 



(-) 



j-\-m— 1 



j — m 



j+m 



2j(2j + 1 
(i - m)(j - m- 1) 



(2j-l)2j(2j + l) 



(C.69) 



(C.70) 



as from which, using the symmetry properties of the 3jm symbols (cf. eq. ()C.55|) . eq. (jC.56l) 
and eq. (|C.57[) ). several other may easily be inferred, for instance 




(C.71) 



1 j-1 
-1 -m+1 



3 

m 



_ r_\-j—m 



(j + to) (j + to — 1) 
(2j - l)2j(2j + 1) 



A number of recursion relations further easing the computation of the 3jm symbols may be 
obtained from the expression displayed in the eq. (|C.82|) of the expansion over 3jm symbols 
of the product of a 3jm symbol with a 6j symbol by giving special values to the arguments 
ii, I2 and l 3 of the 6j symbol and taking help, for instance, from the eq. (|C.88|) to compute 
the involved particular 6j symbol. An illustrating case is the one with the values l\ = \, 
h — js — 5 and Is — 32 — \- Mi can take only the two values ±i and, to comply with the 
zero-sum condition for the 3jm symbols, fi2 must be fixed to /ii +TO3 and /13 to /ii — mi- Thus 



3i 
I is 



32 33 



3i 



33 



31 32 33 
mi TO2 TO 3 



32 



TOi fil + TO3 — fix + TO 2 



/il=±i 



32 



~Ml m 2 Ml - m 2 



2 X 



1 

2 
/'l 



J3 



.73 



-^il - TO 3 TO 3 



(C.72) 



The 6j symbol to compute is the last among those displayed in the eq. (|C.89[) while four among 
the six involved 3jm symbols may be derived, using the symmetry properties (cf. eq. (jC.55l) . 
cq. (|C.56[) and eq. (|C.57|) ). from the one displayed in the eq. ()C69|) . After substitution, we get 



31 

mi 



32 
TO 2 



33 

m 3 



[{32 +m 2 ){j 3 - m 3 )Y 



[{32 +m 2 )(i3 - m 3 )Y- 



3i 

TOi 
31 

m 1 



32 - 

TO 2 

32 - 
TO2 



_L 

2 
1 
2 

_L 

2 



33 ' 

m 3 

33 ~ 
m 3 



_L 

2 
1 
2 

_L 

2 



(C.73) 
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where j — j% + j 2 + 33. Similarly, by setting \\ = 1, Z 2 = js — 1 and I3 = j 2 , it may be shown 
that 



to" + 1)0" -%)0'-%)(7 "% + !)]*( ^ i i) = 
[(j 2 - m 2 )(j 2 + m 2 + l)(j 3 + W3)(j3 + m 3 - 1)] 2 

- 2m 2 [(j 3 + m 3 )(j3 " ma)] 1 f Jl 32 j3 ~ M 
lw /w " V mi m 2 TO3 y 

- [0a + m 2 )(i 2 - m 2 + l)(j 3 - m 3 )(j 3 - m 3 - 1)]* ( ^ ^ x * j (C.74) 

Use can be made of this recurrence relation to compute, for instance, the frequently occurring 
3jm symbol with all the projection quantum numbers m null. If 3 = j\ + j 2 + 33 is odd then by 
the symmetry described in the eq. (|C57j) ) this 3jm symbol is null. If j = ji + 32 + 33 is even 
then from the recurrence relation displayed in the eq. (|C.74[) . using the symmetry specified in 
the eq. (|C.56|) and the eq. (|C.57|) . we may write 



( h h 33 - 1 \ 
I mi m 2 + 1 m 3 — 1 J 



Jl 32 J3 





J2(j2 + l)j 3 (j3 - 1) 



(j + l)(j - 2j 1 )(j - 2j 2 )(j - 2j 3 + 1) 
C7 - 2j 2 - 1)0" - 2js + 2) 



(j - 2j 2 )(j - 2j 3 + 1) 



jl 32 33 - 1 
1 -1 



jl 32 + 1 33-I 




(C.75) 



and, iterating this k times, 



Jl 32 33 





U-2j 2 -2k)\ 



2 fc (^-32)! 



f 2 fc (4-j 3 +fc)! \ 

I Cf-J3-*)I J 



l(j-J2-fc)!j 

(j-2j 2 )!(j-2j 3 )! 
(j-2j 2 -2fc)!(j-2j 3 + 2fc)! 



(j-2j 3 +2k)\ 



u-2] 3 y. 



2 k (^-j 3 + k) 
(i-33- k )< 



ji 32 + k j 3 -k 




$ - h - k)\(i - h + k) \ ( ji j 2 + k 33 - k 

a-32m-33)\ 



whence, setting k = j 3 — i and making use of the eq. (|C.67I) , 



3i 32 3s 




= R i/2 



(3- 2ji)!(j - 2j a )!(j - 2j 3 )! 



(7 + 1)1 



(I)! 



(C.76) 



(|-7i)!(|-i 2 )!(|-i 3 )! 



C.4 6j Symbols and 9j Symbols 

The addition of n > 3 angular momenta Ji + J 2 + • • • + 3 n = J is conceptually not complex, 
but now there is not a single way to build up the coupled states. Let us consider the case 
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n = 3. We can either couple the states \J\M\) and |J 2 M 2 ) to give the resultant \J12M12) then 
couple this with \ J 3 M 3 ) to finally give \ JM) (we can drop the 77. without loss of generality) : 

\{JiJi)Ji2^JM) = J2 \(JiJ2)Ji2M 12 J 3 M 3 ) (JuM 12 J 3 M 3 \J 12 J 3 JM) = 

Mi 2 ,M 3 

= J2 \JiMiJ 2 M 2 J 3 M 3 ) (Ji 2 M 12 J 3 M 3 \J 12 J 3 JM) (,hM x J 2 M 2 \,hJ 2 J 12 M 12 ) (C.77) 

Mi 2 ,M 3 
Mi,M 2 

or alternatively couple the states | J2M2) and | J 3 M 3 ) to give the resultant | J 23 M 23 ) then couple 
this with I J\M{) to again give | JM) : 

\Ji(J2J 3 )J2 3 JM) = \JiM 1 {J 2 J 3 )J2 3 M 23 ) {JiNhJ 23 M 23 \JiJ 23 JM) = 

Mi,M 23 

= ^ |JiM 1 J 2 M 2 J 3 M3<JiM 1 J 23 M 23 |J 1 J 23 JM) (J 2 M 2 J 3 M 3 \J 2 J 3 J 23 M 23 ) (C.78) 

Mi ,M 23 
M 2 ,M 3 

The transformation which connects the states in these two coupling schemes can be written as 
\Ji{J 2 h)J2sJM) = \(JiJ2)Ji2J 3 JM) {{JiJ 2 )Ji2J 3 JM\.h{J 2 J 3 )J 23 JM) (C.79) 

Operating on both side of this equality with J + i it is found out that the transformation 
coefficients ((JiJ 2 )Ji 2 J 3 JM\Ji(J 2 J 3 )J 23 JAl) are independent of M. It then is customary to 
omit it and write these in the form 

v h+j 2 +j 3 +J ) h J J12 



((JlJ 2 )Jl 2 J 3 J|Jl(J2J3)J23J) = [Jl2] [J 23 ] (-) Jl + J2 + J3 + J j j j (C.80) 

L J\ <H J23 J 

which defines the 6j symbol. Inserting the equations of the two coupling schemes into the 
transformation equation which connects them, equating the coefficients of | ,J\M\ J 2 M 2 J 3 M 3 ) 
and using the 3jm symbols in place of the Clebsch-Gordan coefficients, it follows that 



y2(-) j2 ~ j3+M3s+Ji ~ j2s+M [J23] [J] 1 J "' '''' J 1 1 J ' />,; '' 



M 2 



2 -J3 <J23 \ / <J\ J23 

M 2 M 3 -M 23 I \ M x M 23 -M 



1 jl j2 j23 

v I \.h-.h + M 12 + J 12 -.h + M \ j ] \ j] f Jl J 2 ^12 \ f Jl2 J3 J 

X ( } [Jl21 [Jl { M x M 2 -M n J { M 12 M 3 -M 

Using the symmetry properties of the 3jm symbols (cf. eqs. (|C.55j) - (lC.57j) ) this, after appro- 
priate notation substitutions, is more conveniently rewritten 



EH' 1 



+ ' 2 +Ml+M2 



ji k 33 \ f h h 33 
mi m 2 m 3 J \ /ii -/i 2 m 3 



yH^Kilf f f h h h ) ji h h ) (c.81) 

^ K ' NJ 1 'i h h J V -Mi m 2 fi 3 J \ mi -/x 2 -fi 3 J y ' 
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Multiplying both sides by (— ) Ml2 ( "II ^ | and summing over Mi and M 2 , 

\ All AA% — IVl-i 



J 3 
Ji J 2 



J 7' 

J 12 



'23 



7' 

J 12 


h 


J 


M[ 2 


M 3 


-M 


J3 


Ji 


J23 


M 3 


M 2 


-M 23 



<12 

E 

Af 23 ,Mi,M 2 



(-) 



J l + J 2 + J 23 -M l +M 2 +M 23 



Jl 

Mi 



J J23 
-M M 23 



Jl J2 J'12 

-Mi -M 2 M{ 2 



by making use of the orthonormal conditions and symmetry properties of the 3jm symbols. 
After appropriate notation substitutions, this rewrites 



Ji J2 J3 
^1 h h 



Ji J2 J3 

TOl 777 2 771 3 

jl ^2 

7771 M2 



= E h' 1+ ' 2 

Ml, 7*2,7*3 

'3 \ ( h h h 



+(3+/il+A«2+7 t 3 



/l 
/'l 



^2 33 

-^2 m 3 



(C.82) 



It immediately is inferred that the 6j symbol is null unless the four triple of angular momenta 
(jl, 32,33), {31J2 J3), (h,h,h) and (I1J2J3) all satisfy the triangular condition. Inserting 
in the eq. (|C.82[) the algebraic form of the Clebsch-Gordan coefficients gives a complicated 
expression even though to determine the 6j symbol we can set m\ = j\ and 7772 = —ji. 
According to [18] this should give 



Ji J2 J3 
h h h 



E 



A(j 1 ,3 2 ,j 3 )A(j 1 ,l 2 ,l 3 )A(l 1 ,j 2 ,l 3 )A(l 1 ,l 2 ,j 3 ) x 

{-Y(z + i)\ 



(2 ; -jl-j2-j3)K z -Jl- ( 2-i3)!(z-il-j2-i3)!(^-il-i2-i3)Kjl+j2+il+i2-2)Kj2+j3+i2+(3-z)!x 

x(j3+ji+h+h-z)\ 



where 



A(a, b, c) 



(a + b - c)!(a - b + c)\(b + c- a)\ 



{a + b + c+l)\ 

An advantage of this general expression is to put into evidence the invariance of the 6j symbol 
by any permutation of its columns: 



Ji J2 J3 
h h h 



J3 Ji J2 
h h h 



J2 Ji J3 
h h h 



(C.83) 



and by interchange of any of the upper and lower arguments in each of any two of its columns 

h h J3 
ji 32 h 



Ji J2 J3 
h h h 



(C.84) 



It is not be a difficult matter to show that the coupled states of the angular momenta form 
an orthonormal and complete set in the state space of the quantum system under concern and 
that the transformation connecting the states of the different coupling schemes are unitary. It 
then is inferred, using the eq. (|C~80)) ), that 



{{jl32)j' 333-i\(jl32)3" 333i) = ((jlj 2 )j'j3j4 

= E [f\ \f] 



E I j' 1 )jji) ( jl ( J2 33 ) j J4 1 



I(jlj 2 )j"j3j4) 



»1 J Jl J2 J 
33 J4 j 



Jl J2 J 
J3 J4 j 



= S 3',j" 



(C.85) 
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or else from 



{{hh)h2h3\h{jzh)hij) = {{jih)ji2hj\ 



E \jl ih h) J23 j) (jl ( 3233)323 j I 



|j2 0'3jl)j3lj) 



and from 



({3i32)h2h3\Hhji)hi3) = (-y i+32 ^ 12+n+j3 ^ 13 (U23i)3i2333\32(ji33)33ij) 

(Jl(j2j3)j23j|j2(j3jl)j31 j) = (-) jl+j23 ^ {(j2j3)j2S3lj\hUah)331 j) 



(cf. eq. flSSjJ)) that 



E(-) J12+J23+J31 b'23] 



Jl J2 Jl2 
J3 j J23 



J2 33 323 
jl j 331 



33 3l 331 
32 j 312 



(CM) 



Although the general algebraic formula might be used to numerically compute any 6j symbol 
it reveals more convenient to work out simpler expression for special values of the arguments. 
A nice example is provided by setting j'3 = l\ + I2 = UI3 in the eq. (|C.82[) . which then writes 



ji 32 h + h 
h h h 



ji 32 h + h 

TOi TO 2 l\ + I2 



ji h 

TOl fj, 2 



h 
-M3 



- E (-: 

h 32 h 

-/Jl TO 2 /L* 3 



h +h +'3+^1+^2+^3 



h 

1 1.1 



h 

-A*2 



^1+^2 
^1+^2 



(C.87) 



The last 3jm symbol is null unless the zero-sum condition fii — H2 = ~h — h is satisfied, 
but \fii\ < l\ and \^\ < h so that necessarily fii = —li and /J2 = h- On further setting 
toi = — ji the zero-sum condition on the other 3jm symbols impose that TO2 = ji —I1—I2 and 
/j 3 = — ji + I2. The summation in the right hand side then is reduced to a single term and, 
using the symmetry properties of the 3jm symbols, we may write 



ji 32 h + h 



(-) 



h h h 

h h ji 
h ~h + ji -ji 



h+h+h+h ( 3i 



32 



h+h 



= (-) 



-ji-h+h 



ji -ji +h + h -h - h 

h 32 h \ ( h h h + h 

h -h - ji - h ji + h J \ -h -h h + h 



The 3jm symbol in the l.h.s and the two first 3jm symbols in the r.h.s. can be computed using 
the eq. (|C.66[) while the third 3jm symbol in the r.h.s can be computed using the eq. (|C.67|) 
to finally give 



ji 32 h + h 
h h h 



(2ii)!(2i 3 )!(ji+i2+ii+J2+l)Kil-i2+h+i2)!(-ii+J2+ii+i2)! 

(jl-/2+/3)!(j2-/l+/3)! 



(2l 1 +2l 2 + l)\(3 1 +32~h-h)Kji+h-h)K-ji+h+hy-(h+j2-l3)Kh-j2+l 3 y- 
(ji+h+h+iy-ih+h+h+iy- 



(C.88) 
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A few particular instances of this formula are 

= (-y [h 1 ^ 1 ] u =31+32+33) 

20- + i)0'-2ji)0'-2i2)(j'-2i3 + i) 



3i 32 33 

h h 



Ji 32 33 
1 h - 1 32 



3i 



32 



\ 32 



= {-y 

33 



(2j 2 )(2j 2 + l)(2j 2 + 2)(2j 3 - l)2is(2i 3 + 1)_ 
C?-2ja)C7'-2i3 + l) 



32 



3i 

I h - I J'2 



= (-) J 

J3 



(2j 2 + l)(2j 2 + 2)2j3 (2.73 + 1) 

= (_)iT 0' + i)(j'-2j0 



(2j 2 )(2j 2 + l)2j 3 (2j 3 + l) 



(C.89) 



Let us now consider the coupling of n = 4 angular momenta. We can either separately couple 
on one side the states \ J\M\) and | J 2 M 2 ) to give the resultant | Ji 2 Mi 2 ) and on the other side 
the states |J 3 M 3 ) and \J4M4) to give the resultant | J34-M34) then couple the states |Ji 2 Mi 2 ) 
and I J34M34) to finally give | JM) or else separately couple on one side the states \ J\Mi) and 
I J3M3) to give the resultant IJ13M13) and on the other side the states |J 2 M 2 ) and IJ4M4) to 
give the resultant | J24M24) then couple the states | </i3Mi 3 ) and | J24M24) to finally give \ JM). 
The transformation coefficient which connect these two coupling schemes may be written in 
the form 

((JlJ 2 )Jl2 (J 3 J 4 )J34 J I (JlJ 3 )Jl3 (J 2 J 4 )J24 J) = 

( Jl J2 Jl2 ) 

= [J12 ^34 J13 J24] < J3 Ji J34 > (C.90) 

[ Jl3 ^24 J J 

which defines the 9j symbol. It is clear that other coupling schemes can be anticipated giving 
rise to different transformations coefficients and the coupling by pairs of angular momenta we 
have described can also be viewed in terms of three angular momenta couplings, which provides 
with relations between transformation coefficients, for instance 

<(JlJ 2 )Jl2 (J 3 J 4 )J 3 4 J I (JlJ 3 )Jl3 (J 2 J 4 )J24 J) = 

= (( J l J 2)Jl2 (J 3 J4)J 3 4 J I Jl{J 2 (J 3 J4)J34} J234 J) X 

^234 

X (Jl {^2(^3^4)^34 }^234 J I Jl {^3(^2^4)^24 }^234 J) X 

X (Jl {J3(J2Ji)J2i} J234 J I [JlJ3)Jl3 {J2Ji)J2i J) = 

= [Jl2 J34 Jl3 J24] 

p-)*-[4j{* * t}{i t t}{t 1 1} < c - 91 ' 

J 2 34 

(cf. eq. (|C.80p ). Accordingly, after appropriate notation substitutions, 

( hi T22 £ ) = £i r [z 2 ] ( ¥ J ^ J31 ) ( ^ ) ( » ^ ^ ) 

I 331 332 333 J - 1-732 J33 * / \ 321 Z J 23 / \ z Ju Jla / 

(C.92) 
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Using the appropriate closure relations (cf. eq. (IC.40[) ). it may be shown that 
I (Ji J 2 )Ji2 (J 3 J 4 )J34 JM) = I ( J i J 2)Ji 2 M 12 (J 3 J±)J 3i M 3A ) x 

x ((JiJ 2 )Ji 2 M 12 (J 3 J 4 )J M M 34 I (JiJ 2 )Ji2 (J 3 J 4 )J 34 JM) = 

\JiM x J 2 M 2 ) (JiM x J 2 M 2 \{JxJ 2 )Ji 2 Mi 2 ) x 

Mi2,M34,Mi,M 2 ,M 3 ,M4 

IJ3M3J4M4) {J 3 M 3 ,h]VU\(J 3 ,h)J 3i M 34 ) x 

x ((JiJ 2 )Ji 2 M 12 (J 3 J 4 )J 3i M u I (JiJ 2 )Ji2 (J 3 J 4 )J34 JM) (C.93) 

and similarly for | (Ji J 3 )Ji 3 (J 2 Ja)J 2 a JM) then, computing the scalar product of the two 
vectors, using the eq. (|C54|) to replace the Clecbsch-gordan coefficients by 3jm symbols and 
performing the appropriate notation substitutions, 

jn m Ji3 

321 322 J23 
J31 J32 J33 



E 



All 



Jll Jl2 Jl3 \ / J 2 1 J2 2 J23 \ I J 3 1 J32 J33 

TOn mi2 m i3 y ^ m 2 i m 22 m 23 / \ m 3i m 32 m 33 

Jll j 2 l 331 \ ( hi hi 332 \ ( 313 323 333 

mn m 2i m 3 i / I m 12 m 22 m 32 J \ m 13 m 23 m 33 



(C.94) 



A merit of this formula is to immediately reveal the symmetry properties of the 9j symbol. 
This is invariant under even permutations of the columns or of the rows and by swap of the 
lines into columns and the columns into rows. It is multiplied by the phase factor (— )^»,/3 J"./ 3 
under odd permutations of the columns or of the rows. Similarly as for the 6j symbols the 
unitary nature of the coupling schemes lead to the orthogonality relation 



E 

J12 j'34 



Jl 


J 2 


Jl2 ) { Jl 


J 2 


Jl2 


Ji 


Ji 


J34 ) { J3 


Ji 


J34 


Jl3 


J 2 4 


J J I J'l3 


J' 

J 24 


J 



A number of sum rules and recursion relations derived from these can be obtained for the 6j 
symbols and 9j symbols by considering various coupling schemes. Unfortunately, the expres- 
sions often get rather complicated, involving so many parameters that it is difficult to find their 
interrelations and to extract the symmetry of the formulas through the appropriate notation 
substitutions. It gets of course a further puzzling task to analyze the coupling of n = 5, • • • 
angular momenta, which leads to 3(n-l)j symbols. Graphical methods were proposed that 
allows dealing with more ease with the inherent difficulties of this poorly intuitive algebra [17]. 

D SPHERICAL IRREDUCIBLE TENSOR OPERATORS 
D.l Definition and Examples 

A tensor operator T with respect to a symmetry group G is a set of operators that satisfy 
the equivariance property 

VgeG U^TtUig)- 1 = J^T^Cff) (D.l) 



Prl-52 



JOURNAL DE PHYSIQUE IV 



where U(g) <E Q is an image of g in the group Q of the state transformations associated 
with G and T>(g) its matrix representative with respect to a given basis in the state space 
of the quantum system under concern. A spherical irreducible tensor operator is a tensor 
operator that shows the equivariance property for the group of rotations g = R(u,u>) with 
respect to matrix representatives T> J (g) of an irreducible representation T> J of the Lie group 
Q = SU(2). It then consists in a set of 2 J + 1 operators and is naturally labeled by the 
index J. Since any finite rotation can be built from a succession of infinitesimal rotations it is 
equivalent to consider the equivariance property in terms of infinitesimal generators and matrix 
representations of the su(2) Lie algebra. It then may be stated that a set of 2k + 1 operators 
Tq^ (q = —As, — k + 1, ■ • • , k) form the components of a spherical irreducible tensor operator 
T( fe ' of rank k if and only if these operators satisfy the commutation relations 

k 

[J^fHE^WW (m = o,±i) (D.2) 

q' — — k 

more pr ecisely if an only if [ J , T q k) ] = qT q k) and [ J ±1 ,T q k) ] = =p [\{k qp q)(k ±q+l)]* 
(cf. eq. (|C.13p ). It immediately is checked that scalar (k — 0) and vector (k = 1) operators 
arc effectively those defined in the section [C.ll An example of vector operator (k = 1) is the 
differential operator V the spherical components of which, in Cartesian coordinates (as, y, z) 
and spherical coordinates (r,9,ip), are given by 

Vo = d z = cos 9 d r — — — - dg 
r 

V ±1 = T ^(d x ± = jsmfl d r + — do ± — - d v 

Taking the adjoint of the defining relation for T^ k \ calling to mind that U + (a, /3, 7) = 
U _1 (a, ft, 7) and that (ABC) + = C + B + A + , and using the symmetry properties of the 
Wigner D-matrices given in the section [C . 21 it is found that U(a, ft, 7)(Tq fe - ) ) + U~ 1 (a, ft, 7) = 

E 9 '(TfV^;(a^,7) = E 9 '(T? ) ) + (-)«'-«^^_ fl (a J j9 J 7). The set (-)« (T«)+ form the 
components of a spherical irreducible tensor operator (T + )( fe ) called the adjoint to T^ fe ^. 

Let T( fcl ) and Tj( k2 ^ be two spherical irreducible tensor operators wth components and 
Uq^- The quantities 

X <f 5 = E T ^ 1)U S 2) (hqihqzlhfaKQ) (D.3) 
91,92 

transform under rotation similarly as Yl qi q2 \kiqik2q2){kiqik2q2\kik2KQ) = \k1k2KQ) , that 
is according to the irreducible matrix representation D K of SU(2). They thus form the com- 
ponents of a spherical irreducible tensor operator of rank K . It also may be shown more 
directly that the quantities Xq^ do satisfy the required commutation relations. The values of 
K are restricted to K = \ki — fe], \ k% — AS2I + 1, • • ■ , k% + fc 2 and those of Q for a fixed K to 
Q = qi + q 2 = -K,-K+l,--- , K. XW is often denoted {T^ ®V^Y K \ 

The set of {2k x + l)(2fc 2 + 1) operators tJ^U^ 2 ' define the direct product of T( fcl ) and U( fe2 ) 
and can be interpreted as the components of a spherical irreducible tensor operator X^) 



T £0ug») = (kiqik2q2\k 1 k 2 KQ) xg° (D.4) 

K 
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If fei = hi = k then setting K = it is found that 



{T« ® UW} (0) = ^ Tf>U« fc » (fcgfc - g|00) = J2 T ^ )U -^ 



(fc) TT 0) ( k k 
9 -9 



H fe ^H't^u^ = H fe I* -1 ] ( T(fc) ' uW ) ( D - 5 ) 



(cf. eq. ([ClHI) where (T< fe ) • U^) defines the scalar product of and U^. 
If kx = &2 = 1 then setting K = 1 it is found that 

'11 1 

> x;,vu;,y ;j.'7i iu) > t:-,:-'u ;;;"(- )•••• \ >; i 



{ T w (s uw}^ = £t«u« (iftK&iig) = £t«u«(-)V3 



9i 92 



(cf. eq. (|C.54p ). Using the eq. (|C.65|) we compute 



(D.6) 



{t^ ® U« }^ = ^(T.U, - T^U*) = -L(f A U) 2 = -L(f A U) 



V2 V 



where T„ and U m (n, m = x,y,z) are the Cartesian components of the vector operators T 



and U : T x = 

t(i) 



X 1 j — T v 



T (i) T (i) u u; - U {± ' XJ KX) +U KX} 



(D.7) 



U z = Ug 1,1 . It may similarly be computed that 

{tW®UW}^ = -^=(T A U) 

Accordingly, 

(f AU), = i(-)^V6^T gi U 92 ( 1 * J 

91,92 

D.2 Wigner-Eckart Theorem 

Let T( fc ) be a spherical irreducible tensor operator with components Tg . Applying these on 
the angular momentum basis states \tJM), we may build, for each r and J, (2k + 1)(2J + 1) 
states Tg \tJM) (q = —k, — k + 1, • • • , k; M = — J, — J + 1, • • • , J) and from these the states 

i/) T ,j,k(J'M') = J2 T [ k) \ tJM ) (kqJM\J'M') (D.8) 

q,M 

It is an easy matter to invert this, thanks to the orthonormal properties of the Clebsch-Gordan 
coefficients : Y,j'M'( k 9 JM \ J ' M ')( k 9 J1 ^\ J ' M ') = 6 w 6 mm- We § et 

TW |t JM) = Yl A,J,k( J ' M ') {kqJM\J'M') (D.9) 

J'M' 

Using the eq. (|C.13[) and the commutation relations given in the eq. (|D.2p . it may be shown 
that 

3 T^\tJM) = [ Jo,T( fe )]|rJM) + 3 \tJM) = qT^\rJM) + MT^\tJM) 
3 ±1 T^\tJM) = [ J ±1 ,T«]|rJM)+T«J ±1 |rJM) = 



-(fcT#ig + l) 



\tJM) T 



~(JtM)(J±M + 1) 



T^\tJM±1) 
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and, applying J M (/i = 0, ±1) on \J X J 2 JM) = £ MiM2 \J 1 M 1 J 2 M 2 )(J 1 M 1 J 2 M 2 \J 1 J 2 JM), that 

(q + M) (kqJM\j'M') = M' (kqJM\j'M') 



-(k T q)(k±q+l) 



= ± 



(kq±UM\j'M') T 
1 



-(J T M)(J ± M + 1) 



(kqJM±l\j'M'} = 



-(J±M , )(JtM / + 1) 



(kqJM\J'M'T 1) 



Accordingly, 



J±iVv,J,fc(J'M') = T 



JoAMJ'M') = M'i> T! j ik (J'M') 
1 



(D.10) 



■(J'Tl')(/±M'+l) 



A,J,k{J' M ' ± 1) 



which tells that the states ip r .j,k{J' ' M') are either null or proportional to basis states |r' J'M') 
■ ipT,J,k(J' M') — oi t ^j^\t' J'M') where the a Ti j,k are some unknown coefficients that de- 
pends solely on (r, J, k). It then is inferred, since (t"J"M" \ t'J'M') = 5 T " ,t'0~j» ,j'0~M" ,M' 1S 
M" — independent, that 



(t"J"M" I T<*' | tJM) = f{ T ,J,k,T",J") (kqJM\J"M") 



(D.ll) 



This is the essence of the Wigner-Eckart theorem. It allows factorizing the matrix elements 
of the spherical irreducible tensor operators into a geometric component, the Clebsch-Gordan 
coefficient (kqJM\J"M"), and a physical component, the quantity /(r, J, k, t", J"). It actually 
may be derived more elegantly from group theory arguments and is a particular instance of 
the so-called theorem of the three representations. According to this, as applied to the SU (2) 
Lie group, the matrix elements (t"J"M" \ | tJM) under a rotation get transformed 

according to the T> J * <g) T> k (g> T> J matrix representation. The reduction of this cither contains 
the trivial representation once, in which case the matrix element necessarily must not depend 
on the azimuthal quantum numbers, or does not contain the trivial representation, in which 
case the matrix element is null. No other case exists and this is why the geometric component 
factors out. It is customary to express the theorem with a 3jm symbol in the form 



(vjm) f \ T« \(vjm) t ) = (-) 



3f-m s 



jf k ji 
— TUf q TTLi 



{ivj) f || T« || {vj)i) (D.12) 



where ((vj)f \\ |j (vj)i) is called a reduced matrix element. This is often computed in an 
obvious way, for instance from 



(s 



-\ s \s = —m s 



{sf II s II 



we find 



(s f || S || Si ) 

and, more generally, from 

(lm = l\ 1q \lm = I) =1 = 



(D.13) 



I 1 I 

-I / 



(If II 1 II h) 
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we get 

(If || 1 || k) = y/l(l + l)(2l + l) (D.14) 
while from the eq. (|C.63[) we deduce 

(lf\\Y L \\l i ) = ^(-yf{(2lf + l)(2L + l)(l i + l)}^( 1 ^ L Q U A (D.15) 

The Wigner-Eckart shows its full strength in the computing of the matrix elements of spherical 
irreducible tensor operators. Let T^ fcl ) and U^ 02 ) be two spherical irreducible tensor operators 
that act on part 1 and 2 of a system with angular momenta ji and ] 2 . The matrix element of 
{T( fel > <g> U( fe2) }£f 3 = Xjjf 5 between the states \vjij 2 JM) fti ) may be written 



( (vjih JM) S \ Xf | (vhh JM)i > = 

= (-) Jf - Mf (_% q M t ) (( v jd2J)f\\X {K) \\(vji32J)i) (D.16) 

but since T< fcl ) and are supposed acting on different state spaces, it follows that 

((vj 1 j 2 JM)f\X<g ) \(vj 1 j 2 JM) i )= ]T {J f J l ](-yif-»f+ M f+^-^+ M > x 

mi/ ,m2 / ,mii ,rri2i 

i jxf J2f 3 L \( hi 321 *{OW fcl k2 K n ){-) kl - k2+Q \K]x 



m lf m 2 f -Mf J \ mi; m 2l -Mi J ^ \ qi q 2 -Q 



11,12 



x {ivjimxj 2 m 2 )f\ T ( ^U^\(vj imi j 2 m 2 ) f t 

— [JfJiK] (-.yif-32f+ M f+3^~]2i+M i +k 1 -k 2 +Q 

mi f ,rri2f ,m>2i ,11 ,12 ,P 

x f hf hf Jf \ ( jii hi Ji \f h k 2 K 

V m if m 2f -M f J \ mu m 2i -Mi J \ qi q 2 -Q 

xR— /(_^ I ^ ) ((^HT^IIp^) x 
xR— £ * ) (p^lluWlK^),) (D.17) 

On multiplying by (— } J f~Mf ^ ^ ^ ^ and summing over M/, Af^ and Q it may 
be shown that 

({vhhJ)f\\X- {K) \\{vhhJ)) = 

!jif ju ki "1 

hf hi k 2 \ (D.18) 

Jf Ji K J 
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Setting k 2 = the reduced matrix element for the operator T( fcl ) acting only on part 1 of the 
system is given by 

= S j2f , j2i (-) jlf+hf+Ji+kl [JfJi}{ J / - 1 -\ \{{vh)f\\^ ) \\{vh)) (D.19) 

^ Jli J2f Jlf ) \ ' 

Setting ki = the reduced matrix element for the operator U^ fc2 ^ acting only on part 2 of the 
system is given by 



-S nf , ni ^y i ^ +Jf+k2 [JfJ l ]{2 l f lf }((^)/l|U (fe2) IIN 2 ),) (D.20) 

If there is no possibility to regard T( fcl ) and U^ fe2 ^ as acting on separate parts of a system then 
there is no point in indicating that J is the resultant of ji and j 2 . In this case 

{{vJM) f \X%\vJM) l ) = 

-k 2 +Q tif] ( h k 2 K 



= [tf] f £ 2 ^ U(vJM) f \T^)V^\{vJM) i ) = 

= E _(-p- k ^ [K] ( £ ^ K Q ) ((vJM) f \T^\JpM)(JpM\V^\(vJM) i ) = 

qi,q2,J,p.M 

' ki k 2 K 



_ ^_jk 1 -k 2 +Q+jf-M f +j-M ^ 

91, 92, J, M 
Jf k\ J \ ( J k 2 Ji 

-Mf q x M ) \ —~M q 2 M f 



(ivJ) f \\T^\\pj) (pJ\\V^\\{v.J)) = 

E ( _ ]Jf+ J i+ K+J f -M f , K] f Jf K J t \ ( Jf K J t 1 
{ ' 1 1 \ -M f Q M t )\k 2 J fci J 



((«J) / ||T^)||pj)(pJ||U(*»)||(t;J) i ){ J * j M (D.21) 



so that 



((vj) f \\x^\\pj) = (-yf^[K]j2{ 



J f K J 4 
k 2 J k\ 



x ((t^UT^Upj) [pJ\\V (k ^\\(vJ) t ) (D.22) 



E MANY EQUIVALENT ELECTRON STATES 

It may be shown that an TV— electron states \sf,i) = \{vLMlSMs) /,») of the configuration 
(nZ)jv can be built up from (A/ - — 1)— electron states |<T/^) = \(vLMlSMs) of the configu- 
ration (nl)ss-i and 1— electron state !»?/,») = |(nZm;sm s )/ ; j) of the configuration (nZ) as 

\(vLM L SMs)f,i) = ]T 1} n/a)^) 
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where 

\(vLSnls) u )= Yl E ((^lm)fA(LM L ) fti )((SM^sm s ) f , i \(SM s )f, i )x 

(ML,Ws)f,i 

x \ (vLM L SM s ) f ,i) \{nlmism s )f ti ) 

are the coupled states of (Af— 1) electrons with the Af th electron and where ((vLS) f y i\}(vLS) /^) 
are called coefficients of fractional parentage. This concept is extremely useful when deal- 
ing with matrix elements of sums O = ^2 u O(u) of one-electron operators 0(u). Indeed, 
owing to the equivalence of the electrons, the matrix element of 0(u) between the states 
\(vLMLSMs)f,i) does not depend on which electron u has been chosen to compute it, so that 

((vLM L SM s )f\ O \(vLM L SM s )i) =M((vLM L SM s ) f \ 0(AT) \(vLM L SM s )i) (E.l) 

that is 

((vLM L SM s )f\ O \(vLM L SM s )i) = AfY((vLS) f {\vL^)(vLS\}(vLS)i) x 

vTS 

x {<J^lmi) f \{LM L ) f ){<JMllm l ) i \{LM L ) i )x 

M Lf M L M Li 

x Y ((Wsm s ) f \(SMs)f)((LM^sm s ) i \(SM s )i) x 

MsfJJsMsi 

x ((lmism s )i\ O(TV) {lm l sm s ) l ) = (E.2) 

or else 

((vLM L SM s )f\ O \(vLM L SM s )i) = 

= MY,( f{fi)@\} i)((i SLM L SM s)f \ OM \(0LM L SM s )i) (E.3) 

6 

by introducing self-explicit shorthand notations to simplify the equations. 
If the operator O is a tensor operator of rank k independent of the spin 

TW^^t^fn) (E.4) 

u 

then applying the Wigner-Eckart theorem 
((vLM L SM s )f\ Tf {(vLMlSMs),) = 

- r l - m - ( X f X m Li ) ^ lsm ^\ T(fc) = 

6 

x (j9LSM)f\\ t^(A0 \\(9LSM)ij (E.5) 
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Use can be made of the eq. (|D.19[) to compute the reduced matrix element assuming that the 
part 1 of the system is played by the Af th electron. We get 

((vLM L SM s ) f \ T« {(vLMlSMs).,) = 



The operator T( fel,fe2 ) comprising the (2fci + l)(2fc 2 + 1) components Tq^^ is a double tensor 
if it behaves as a tensor of rank ki with respect to the spin angular momentum states and as 
a tensor of rank k% with respect to the orbital angular momentum states. If this is the sum of 
one-electron double tensors, spin (x) - orbit (y) direct products, then 

u u 

Its matrix elements are given by 

((vLM L SM s ) f \ \(vLM L SM s )i) = 

_ (_\Sj-M S f+Lf-M L} f S f k x Si \ ( Lf k 2 £i \ 
1 1 \ -M Sf qi M Si ) \ -M Lf q 2 M u ) 

x AfY / (0f{\O)(O\}O l ) ((8LSM) f \\ t^' k ^{Af) \\(§LSM)i) (E.8) 



Use again can be made of the eq. (|D.19j) to compute the reduced matrix element applying it 
separately to the spin and to the orbital part of the double tensor with the part 1 of the system 
played by the M th electron. We get 

{{vLM L SM s )}\ T^'^ \(vLM L SM s )i) = 

L f k 2 Li 



(_\Sf-M S f+Lf-M L f ( S f k x Si \ / 

1 ; V - M sf qi M s l ) \ 



-M Lf q 2 M Ll 

X 



where 



: AA^(0 / {|0)(6?|}^)(-) S+L+s+ ' +S+L+fcl+fc2 [SfSiLfLi 

(sl\\ t (fel ' fc2) (AA) \\sl) (E.9) 

(sl\\ t^ M \Af) \\slj = («|| x (kl) (Af) \\s) y (fc2) (A0 ||A (E.IO) 



Sf h Si \ ( Lf Li 
s S s | 1 I LI 
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